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ABSTRACT 

We compare atomic gas, molecular gas, and the recent star formation rate (SFR) inferred from Ha 
in the Small Magellanic Cloud (SMC). By using infrared dust emission and local dust-to-gas ratios, we 
construct a map of molecular gas that is independent of CO emission. This allows us to disentangle 
conversion factor effects from the impact of metallicity on the formation and star formation efficiency 
of molecular gas. On scales of 200 pc to 1 kpc (where the distributions of H2 and star formation match 
well) we find a characteristic molecular gas depletion time of r"^°' ^^ 1.6 Gyr, similar to that observed 
in the molecule-rich parts of large spiral galaxies on similar spatial scales. This depletion time shortens 
on much larger scales to ~ 0.6 Gyr because of the presence of a diffuse Ha component, and lengthens 
on much smaller scales to ^ 7.5 Gyr because the Ha and H2 distributions differ in detail. We estimate 
the systematic uncertainties in our dust-based r"^°' measurement to be a factor of ~ 2-3. We suggest 
that the impact of metallicity on the physics of star formation in molecular gas has at most this 
magnitude, rather than the factor of ^ 40 suggested by the ratio of SFR to CO emission. The relation 
between SFR and neutral (H2 -I- H I) gas surface density is steep, with a power-law index k, 2.2 ± 0.1, 
similar to that observed in the outer disks of large spiral galaxies. At a fixed total gas surface density 
the SMC has a 5 — 10 times lower molecular gas fractio n (and star formation ra te) th an large spiral 
galaxi es. We explore the ability of the recent models by iKrumholz et al.l (|2009c[ ) and lOstriker et al.l 
()2010t ) to reproduce our observations. We find that to explain our data at all spatial scales requires 
a low fraction of cold, gravitationally-bound gas in the SMC. We explore a combined model that 
incorporates both large scale thermal and dynamical equilibrium and cloud-scale photodissociation 
region structure and find that it reproduces our data well, as well as predicting a fraction of cold 
atomic gas very similar to that observed in the SMC. 
Subject headings: galaxies: ISM — ISM: clouds — galaxies: dwarf, evolution — Magellanic Clouds 



1. INTRODUCTION 

The relation between gas content and star formation 
activity in galaxies has been a matter of intense inves- 
tigation since the pioneering work of Maarten Schmidt. 
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iSchmidtl (|1959( ) suggested that the star formation rate 
(SFR) in a galaxy is proportional to a power of the gas 
density, such that /Osfr '^ Pgasj where n ~ 1 — 3 and 
most likely n ~ 2 based on a number of arguments that 
included the luminosities of open clusters, the abundance 
of helium, and the vertical distribution of objects in the 
plane of the Milky Way. Here we refer to the quantita- 
tive relation between gas and star formation as the "star 
formation law" for convenience without intending to sug- 
gest a physical law or a specific functional form. 

Most modern empirical studie s of the ext r agalactic sta r 
formation law follow those bv iKennicuttI (|1989l . |1998[) . 
who studied disk-averaged correlations in a sample of ~ 
100 galaxies including starbursts and high-mass dwarfs. 
This work linked the surface density of star formation 
rate, SsfRj to the surface density of tot al neutral (atomic 
and molecular) gas, Egas = Shi + Sh2- iKennicutg (|1998[ ) 
found that Esfr oc Sgas^^^ with I + p k. 1.4 ± 0.15 for 
his composite sample of galaxies. Since the gas deple- 
tion time rl^p = Sgas/SsFR k Egas~^ for a power-law 
relation, p = indicates that the gas depletion time is 
constant (independent of environment), while p > in- 
dicates that star formation is more rapid in high-Egas 
regions. In practice, it is also important to keep in mind 
the observational complexities associated with measuring 
Egas and particularly Eh2 from CO observations, and the 



systematics thus introduced. 

Several subsequent studies focused on the star for- 
mation law within galaxies, employing high resolution 
molecular and H I observations ([M artin fc Kennicuttl 
dOOllWong & Blitz"2002- Boissicr ct al. 2003). More re- 
cently, Bigicl ct al. (200 8) and Lcrov ct al. (2008 ) ana- 
lyzed 12 nearby spirals with high quality H I, CO, far- 
infrared, and far-ultraviolet data. These observations 
show a clear correlation between Ssfr and the sur- 
face density of molecular gas, Sh2, with an approxi- 
mately constant star formation rate per unit molecular 
gas, yielding an approximately linear power-law index 
l-t-p = 1.0±0.2. By contrast, they found a steep cor- 
relation between Ssfr and the surface density of atomic 
gas, Shi (the observed correlation in the optical disks of 
large galaxies has 1 + p > 2. 7 — 3.5, and 1 + p ^ 2 i n 
the H I-domin a ted ou ter disks Bigiel et al. 2008, ■2010b). 
ISchruba et"aLl (j2011[ ) and Bigicl ct al. (2011) extended 
these findings to a large r sample of 30 d i sk gal axies, while 
iBlanc et alT(|2009D and iRahman et~aLl (|2011[ ) arrived at 
similar conclusions in very detailed studies of individual 
targets. 

A consequence of the observed linearity of the star 
formation law over the molecule-dominated regions of 
disks is that the typical time scale to deplete the molec- 
ular gas by star formation in a disk galaxy is t^°^ = 

Smoi/SsFR ~ 1.9 ± 0.9 Gyr (Emoi corresponds to I]h2 
corrected by a 1.36 factor to account for the cosmic abun- 
dance of helium) . The lack of dependence of t^"^ on en- 
vironment can be naturally understood if two conditions 
are fulfilled. The first condition is for star formation to 
be a local process primarily determined by the condi- 
tions inside the giant molecular clouds (GMCs) where 
it takes place. GMCs are themselves isolated from the 
global galactic environment provided that they are self- 
gravitating and ther efore overpressured with respect to 
the neighboring gas ()McKee fcOstrikeri [20071 ). The sec- 
ond condition is for the properties of GMCs to be uni- 
versal, and therefore independent of their galactic en- 
vironment. Because star formation takes place in the 
densest regions of GMCs, themselves self-gravitating and 
thus mostly decoupled fr om their surroun dings, the first 
condition appears likely ([Krumholz fc Mc Kcc 2005), at 
least in mid-disk and outer-galaxy regions. In galactic 
centers it is less clear whether there are isolated, gravita- 
tionally bound GMCs, or instead a distributed molecular 
medium. Environmental considerations may be more im- 
portant if cloud-cloud collisions play an important role 
(jTanll200"o( ) . Similarly, the second condition appears to be 
broadly satisfied as the resolved properties of GMCs (or 
at least of the dense regions of GMCs, bright in CO) are 
observed to be remarkably constant over a wide range of 
galaxy environments (the Local Group spirals and nearby 
dwarfs, Bolatto et al. 2008; the Milky Way, Heyer et al. 
2009; the LMC, Hughes et al. 2010; the outer disk of 
M33, Bigiel et al. 2010). 

Most of the studies of the relation between gas and star 
formation, and the ensuing conclusions, focus on large 
galaxies which tend to be rich in molecules. There is 
a dearth of comparable information for low-mass low- 
metallicity star-forming dwarf galaxies, primarily be- 
cause such objects emit only very faintly in molecular 
gas tracers such as CO. The measurements that do exist 



reveal large ratios of S FR to CO emission in late-type, 
low-mass galaxies (e.g.. lYoung et al. 1996; Ca rdan et al.l 
[20071 iLerov eFani2007bl:IKrumholz et al.ii201lTr ^ 

Studies of the star formation law in dwarf galaxies are 
interesting for a number of reasons. Fundamentally they 
probe a very different physical regime from large spi- 
ral galaxies, one where atomic gas dominates the sur- 
face density of the interstellar medium (ISM) on large 
scales, and elements heavier than helium are less abun- 
dant. This dearth of heavy elements causes a num- 
ber of differences in gas chemistry and physical condi- 
tions. For example, dust-to-gas ratios are lower, pro- 
ducing lower extinctions and higher photodissociation 
rates (the hkely cause of their general CO faintness, 
Israel et al .lll986HMalonev fc Blacklll988HLeaueux et alj 
1994; Bola tto et al.lll999n . 

More importantly, studies of the star formation law 
in low-mass low-metallicity star-forming dwarf galaxies 
provide fertile testing ground for star formation theo- 
ries. A natural implication of the approximately constant 
molecular tJ^°' among normal galaxies is that the effec- 
tiveness at forming molecular gas plays a critical role 
in establishing the star formation rate. The molecular 
gas fracti on varies systematically both within and among 
galaxies ([Y oung fc ScoviUg [19911: iWong fc Blitj l2002t 
iBlitz fc Ro solowsk\ n2006l : ILerov et al.l l2008f ). Dwarf 
galaxies appear to have low molecular fractions, despite 
their large gas fractions and long time scales to deplete 
gas reservoirs at current rates of star formation. If their 
r^"°' values are similar to normal galaxies, it would sug- 
gest that CMC formation (and destruction) is the rate- 
limiting step for star formation in dwarf systems. 

Because their low metallicities contrast with those in 
spiral galaxies, star-forming dwarf galaxies can break a 
number of degeneracies in the physical drivers of the 
molecular fraction. Different models predict very differ- 
ent behaviors for the impact of metallicity on molecular 
fractions and star formation rates according to their em- 
phasis on dynamics, thermodynamic equilibrium, or dust 
shielding. For example, the requirement o f a minimum 
extinction for sta r formation to occur (e.g., lMcKeelll989l : 
iLada et"aL[|2009D would have a proportionally more im- 
portant impact in low metallicity objects, while models 
that are solely based on dynamics cannot distinguish be- 
tween low and high metallicity. 

The main obstacle to using dwarf galaxies to test var- 
ious aspects of the star formation law is the difficulty 
in observing their molecular gas distribution, stemming 
from the faintness of their CO emission and the uncer- 
tainty in its quantitative relation to H2 . The Small Mag- 
ellanic Cloud (SMC) offers a unique opportunity in this 
regard. Its proximity allows us to probe small spatial 
scales with observations of modest angular resolution. 
As a result, we can use dust emission to construct a map 
of EH2 that does not depend on CO. In this study we 
use such a map to compare the distributions of molecu- 
lar gas, atomic gas, and recent star formation traced by 
ionizing photon production. 

1.1. The Small Magellanic Cloud as a Galaxy 

By virtue of its location and properties the SMC is 
a prime laboratory for the study of the relation be- 
tween gas and star formation in low-mass galaxies. In- 



deed, the SMC was the target of the first extragalactic 
study of the star fo rmation law, rela ting H I and stel- 
lar surface densities ('Sandulcak"1969 [). Locate d scarcely 
61 kpc away fe ilditch ct al. 2005; K eller fc W ood 2006; 
ISzewczvk et al.l 120091) . the SMC is the nearest gas- rich 
low metallicity d warf galaxy w i th active star-formation 
(ZsM c ~ ^0/5, IDufoud [l98l IKurt et alJ HOOl IPaeeil 
|2003[) . As such, it provides invaluable insight into the 
physics and chemistry of the ISM in chemically primitive 
star-forming galaxies. Moreover, because of its proximity 
it is possible to carry out studies on the scale of individ- 
ual young stellar objects, which reveal subtle differences 
in the temperature and chemistry of star- forming cores 
(jvan Loon et al.ll2010HOliveira et al.ll2011[ ). 

With an HI mass o f Mm ^ 4.2 x 10^ Mq 
(|Stanimirovic et all 119990 the SMC is rich in atomic 
gas, although als o strikingly defhcient in cold H I 
(jDickev et al.ll2000( ). It was already apparent in the early 
observations that the distribution of the H I is complex, 
with more than one com ponent along the line-of-sight 
(jMcGee fc NewtonI ri981[ ). This complexity is at least 
in part due to the prese nce of several supergiant shells 
(jStanimirovic et al]|1999[ ). The underlying dynamics ap- 
pear to be disk- like with an inclination i w 40° ± 20°, al- 
though disturb ed by the interaction wit h the Milky Way 
and the LMC (jStanimirovic et al]|2004f) . In our analysis 
we will ignore this complexity, but the geometry of the 
SMC and the level of turbulent support in the gas remain 
some of the most important caveats in the comparison 
to star formation models. 

The SMC hosts a healthy amount of star formation 
despite being disproportio nately faint in CO emission 
(jKennicutt et a l. 1995; Isra el etal.lll986lll993D . In fact, 
in the SMC the CO is under-luminous with respect to 
the star formation activity by almost two orders of mag- 
nitude when compared to normal disks, or even more 
massive small galaxies. A possible explanation for this 
fact is that the SMC is extraordinarily efficient at turn- 
ing the available molecular gas into stars (i.e., the II2 
depletion time is very short), which would suggest that 
we are either observing an out-of-equilibrium situation (a 
fleeting starburst) or that the conditions conspire to keep 
a small reservoir of extremely short-lived GMCs. A more 
likely alternative is that the weak CO emission is not rep- 
resentative of star-forming molecular gas, with H2/CO 
considerably larger than in the Milky Way ([Israel et al.l 
[1986; Rubi o ct al. 1993). 

iLerov et al. (2007^ used new far-infrared observations 
to map the H2 distribution in the SMC bypassing CO 
emission, following an extens i on of t he methodology pre- 
viously employed by llsraell ([1997t ) in the Magellanic 
Clouds. The feasibility of using dust to map II2 was 
clearly demonstr ated bylDame et al.l ([2001D in the Milky 
Way (see also iBloemen et al.l 119901 for an early study 
of the correlation between dust emission and gas in the 
Milky Way). It was also sho wn to be c onsistent with 
virial masses on large scales ([Lerov et a l. 2009). Fur- 
thermore, a systematic application of an extension of 
this methodology throughout the Local Group produces 
molecular masses that are consistent with those obtained 
by other methods at approximately Galactic met a llicities 
([Lerov et al.ll20lH ). The analysis bv ILerov et all ([2007a[ l 
shows that the column density of molecular gas present 



in the SMC is far in excess of that derived by applying 
the Galactic CO-to-H2 conversion factor, Xqo: to the 
CO observations. Their modeling yields a total molec- 
ular mass Mh2 ~ 3.2 x 10^ Mq. Therefore the neu- 
tral ISM in the SMC is approximately 10% molecular, a 
lower fraction than o bserved in normal late spirals but 
not dramatically so (lYoung fc Scovilldll991l). Fur ther 
refinements to the analysis ( Lerov et al.ll2009l 120111 and 



this work) broadly confirm these numbers. 

In this study we analyze the spatially resolved correla- 
tion between the distributions of molecular gas, atomic 
gas, and recent star formation traced by ionizing photon 
production on different scales, comparing the star forma- 
tion law in this small low-metallicity galaxy with that in 
large disks. We present our data and discuss the method- 
ology we use to measure molecular column densities and 
star formation rates in 321 We show our main observa- 
tional results in SjSj focusing on the relation between II2 
and star formation in H3.21 total gas and star formation 
in ^^ In 21 we compare our results to recent analytical 
physical models of star formation in galaxies. Finally, we 
summarize our conclusions in ^Sl 

2. OBSERVATIONS AND METHODOLOGY 
2.1. Estimating H2 from Infrared Dust Emission 

To estimate molecular gas surface densities from far- 
infrared emission, we use an extension of the methodol- 
ogy already employed by ILerov et"aLl ([2009D . Here we 
summarize the main points, but the reader should refer 
to that paper for the details, as well as a more thorough 
discussion of the effects of several systematics involved 
in the derivation of the H2 map. 

We combine H I and infrared (IR) imaging to esti- 
mate the distribution of H2 at a resolution ~ 40" 
(r ~ 12 pc). The IR data originate from the combina- 
tion of two Spitz er surveys, SAGE-S MC ([Gordon et al.l 
[20TI and S^MG (IBolatto et al.l [200l. The H I data are 
from 'Stanimir ovic et al.l ([1999f) and include both inter- 
feronietric (from the Australia Telescope Compact Array, 
ATGA) and single-dish (from Parkes) observations, and 
so are sensitive to emission on all spatial scales. When it 
is necessary to include an inclination correction we use 
i « 40° for the SMC, as derived from fitting the H I kine- 
matics ([Stanimirovic et aT]|2004D . Our analysis uses data 
from Spitzer^s MIPS instrument at 70 and 160 /im. We 
correct the 160 fim map for foreground contaminating 
emission from Galactic cirrus by subtracting the appro- 
priately scaled H I emission from the Milky Way, which 
is identified on the basi s of its velocities (see iBot et ahl 
l2004t [Lerov et aI1l2009f ). Our gas surface densities in- 
clude a 36% correction for the mass contribution of he- 
lium. 

We derive the surface density of H2, Sh2, by model- 
ing the IR emission to infer the amount of dust along 
each line of sight. Together with a dust-to-gas ratio (or 
strictly speaking a dust optical depth to H I -f- H2 surface 
density ratio) , which we estimate by comparing dust and 
H I iteratively (see below), the dust emission allows us 
to estimate the total amount of gas present. Subtract- 
ing the measured H I surface density yields SH2 ■ We 
should caution that this procedure may identify as H2 
very cold self-absorbed H I gas possibly associated with 
cloud envelopes, simply because this material will be dis- 
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Fig. 1. — (Top) H2 column density map at ~ 12 pc resolution. This map is obtained from modeling the Spitzer dust continuum 
observations from S^MC/SAGE-SMC IIBolatto et al.|[2n0^ : [Gordon et al.ll201in together with the combined ATCA/Parkes 21 cm H I map 
l|!Jtanimirovic et al.lll999l ). The colorbar inset indicates the values for the color scale, in units of 10^^ cm~^. The Nh2 contours are placed 
at Njj2 ^ 1.4, 3.4, 8, and 12 X 10^^ cm~^, equivalent to deprojected molecular surface densities Sj^ol ~ 23, 56, 130, and 200 M0 pc~^ when 
including the correction for the cosmic abundance of helium and the 40° inclination of the source. The western region, oriented roughly 
SW to NE and harboring most of the star formation activity as well as the molecular gas is called the SMC Bar. The extension to the SE 
is called the SMC Wing, and is unremarkable in molecular gas except for the N83/N84 molecular cloud complex, which is the Eh2 peak 
of the galaxy. (Bottom) H2 map overlaid on the unobscured Ssfr map derived from Ho. The SFR is computed using the first term of 
Eq. |4] and the color scale ranges from Ssfr = to 0.1 Mq yr~^ kpc"'^. There is good correlation between regions with star formation 
and regions with molecular gas, but not a one to one correspondence at high spatial resolution. The overall correlation improves when 
smoothing to scales of 200 pc and larger. There is also a pervasive component of diffuse, low level Hq emission. 



proportionally faint at 21 cm and thus incorrectly sub- 
tracted. We include a correction for optic al depth that 
accou nts for self- absorption in the H I map (|Dickev et alj 
I2OOOD , but such correction is by necessity only statistical. 
Thus the procedure we outline may be more accurately 
thought of as producing a map of the very cold neutral 



and molecular phases, which under normal circumstances 
would be completely dominated by the mass of H2 . 

We use the optical depth at 160 /im, rigo, as our proxy 
for the amount of dust along the line of sight. We es- 
timate this quant ity following the method outlined by 
iLerov et al.l ()2009f ) , though note that our calculations im- 



prove on theirs in several aspects. We take the ratio of 
IR intensities at 100 and 160 /xm {Iioo and /leo) to be 
a tracer of the equilibrium temperature of a large grain 
population that contains most of the dust mass. Because 
the 100 fiia intensity comes from IRAS, we only known 
it at ^ 4' resolution. To overcome this limitation to the 
resolution of the study we compare the /ioo/-^i60 ratio 
to the ratio of IR intensities at 70 and 160 /im. The 
emission at 70 /im, /70, suffers from significant contami- 
nation by very small grains out of equilibrium, and it can 
only be used to determine the temperature of the large 
dust grains after removing the contamination. We cal- 
ibrate a relationship that allows us to predict /ioo/^i60 
from /70//160, at the coarse resolution of the 100 /Ltm 
data. The fact that the iMno/hfio-Iro/h m relation is 
well-defined can be appreciated in Figure 2 in lLerov et al.l 
((2009,i), and the fitted relations (Equations 2 and 3 in the 
same paper) suggest that half of Ijo is due to emission 
from stochastically heated grains. We apply our calibra- 
tion at the full resolution of the 160 /xm data (~ 40") to 
estimate the dust temperature and optical depth, Tigo- 
Thus this procedure allows us to take full advantage of 
the 160 fim (~ 12 pc) resolution while guaranteeing that 
we match the results derived from the 100 and 160 fiui 
on size scales larger than 70 pc. 

We now have images of Tieo, which linearly traces the 
distribution of dust, and the H I column density from 
21 cm radio observations corrected for optical depth ef- 
fects (' Dickev et al.ll2000l ). To create an image of EH2 
we combine these maps to produce local estimates of our 
gas-to-dust ratio proxy, ^gdr — ^gus/Tieo- Note that 
^GDR is an observable in regions dominated by atomic 
gas. We use it because the actual value of the mass emis- 
sivity coefficient for dust is unimportant in producing the 
H2 map. Its systematic large-scale variations, however, 
matter and are very important to consider as we discuss 
below. Hence 



^H2 



Tieo S, 



GDR 



■^HI- 



(1) 



There are several ways to obtain the (5gdr, from maps 
of dust and H I. One could simply take the average ra- 
tio across the entire SMC, which is certainly H I domi- 
nated on large scales. Unfortunately, the SMC is known 
to harbor significant large-scale variations in (5gdr, with 
higher values in diffuse, low densit y regions, and par- 
ticularly in the SMC Winp; (Stanimirovic et al.l 120001: 

Similar large scale 



IBot et all [20041: iLerov et al 



^ i stai n 
[2007al) . 



variations, for example, appear to be present in the frac- 
tion of dust in aromati c polycyclic hydrocarbons (PAHs, 
iSandstrom et al.l l2010^ and in the optical extinction per 
gas column density (Dobashi et al. 2009) between lines- 
of-sight with dense and diffuse gas in the SMC. These 
variations are probably not rooted in variations in metal- 
licity. Indeed, metallicity appe ars to be app roximately 
uniform across the source (e.g.. iDufouH 11984 ). although 
note that these measurements are limited to the dense 
gas phases that form stars. Rather, they may be due 
to grain processing in the ISM. These variations could 
be analogous to the changes in dust optical depth per 
unit mass driven by temperature and grain structure 
changes observed ove r smaller scales in t h e Milky Way 
jBernard et al.' '1999'; 'Schncc ct al.' '20081 iFlagev et all 
[2009; Planck Collab oration et al., 201111). An example 



of a nearby galaxy which appears to harbor such large- 
scale variations is M 31 ([Nieten et al.l 120061) . Evidence 
suggesting such variations is also observed in dwarf ir- 
regulars, part icularl y those with e xtended H I envelopes 
(iWalter et al.ll2007t ' lDraine et al.ll2007[) . In the case of 
the SMC, (5gdr appears to be systematically different in 
the Wing, t he SW region of t he Bar, and the NE region 
of the Bar (jLerov et al.ll20lH) . 

In the presence of systematic Jgdr variations adopting 
a single (5gdr, although appealing in its simplicity, will 
lead to dramatically overpredicting the amount of II2 in 
regions that have an intrinsically lower ^gdr- A more 
robust approach and considerably more conservative is 
to do a local average determination of ^gdr, either as a 
function of either position or some other environmental 
quantity like IR surface brightness. That is the procedure 
we follow, and the H2 results thus obtained are consis- 
tent with other obser vational constraint s, such as virial 
masses on large scales (jLerov et al.ll2009D and II2 column 
densities determined in absorp tion along diffuse lines- 
of-sight (jTumlinson et al.ll2002D . Nonetheless, (5gdr re- 
mains the largest source of systematic uncertainty in our 
estimation of Sh2- It is possible we are missing a large 
extended molecular component, although we consider it 
unlikely. 

We implement this approach, making an iterative local 
measurement of ^gdr- First, we blank regions of the rigo 
and H I maps where bright CO is detected — these lines 
of sight are likely to include II2, so that Sni/Tieo is not 
a good tracer of the total (5gdr- Next we smooth both 
maps with a kernel of radius R ^ 200 pc and calculate 
(^GDR = Shi/ti6o for each region at this coarser resolu- 
tion. This choice of kernel size represents a compromise 
between, on the one hand, calibrating ^gdr as close as 
possible to the molecular regions to minimize the effect 
of spatial variations, and on the other, avoiding includ- 
ing regions with significant II2 in its calibration. Using 
a kernel with R ~ 450 pc increases the total molecular 
mass of the SMC by ~ 40%, while reducing the kernel 
to i? ^ 100 pc reduces the mass by ~ 50%, thus brack- 
eting the impact of the choice of the smoothing scale. 
The (5gdr niap produced this way has a median value 
^GDR ~ 5.9 X 10^^ cm~^, and a spread of ~ 0.4 dex. Rep- 
resentative values for the N83/N84 region near the end of 
the Wing, the SW of the Bar, and the NE of the Bar are 
(^GDR ^ 5.9, 4.8, and 3.1 x 10^^ cm~^ respectively with 
the highest values occuring near the central portions of 
the Wing. These values span the range of 8 to 14 times 
higher than Galactic, taking ^gdr ~ 4.1 x 10^'' cm~^ 
determined for Galactic cirrus as representative of the 
Milky Way (jBoulanger et al.lll996D . 

We apply this low-resolution ^gdr to the nearby 
blanked regions with bright CO emission. We then apply 
Equation [1] to our full-resolution maps to estimate SH2 
everywhere. Because this may reveal new regions where 
H2 makes a significant contribution to the gas surface 
density we iterate the process once, blanking everywhere 
that has bright CO and everywhere with I]h2 > 0.5SHI- 
We verified that iterating further does not significantly 
change the molecular surface densities. The result of this 
second iteration is our estimate for I]h2 at an angular res- 
olution 9 ^ 40" (Figure [J). The map shows very good 
correspondence with the recent optical reddening maps 
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from iHaschkc et al.. (,2011,) . 
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We estimate the uncertainties in this map via a Monte 
Carlo calculation. In each iteration of this calculation, 
we add realistic noise to the data, adjust the zero point 
of the IR maps within the uncertainties, rederive the re- 
lation used to predict /yo/^ieo from /ioo/-^i60 ^-nd propa- 
gate the noise in this relation forward. We also vary the 
wavelength dependence of the dust opacity, /3, across its 
likely range /3 w 1 — 2, with one value randomly chosen 
for each iteration. We carry out 1000 such exercises and 
calculate the 1,2, and 3cr uncertainties in Tdust and rigo 
from the resulting distributions. These are propagated 
into uncertainties in S]h2- Using this technique we esti- 
mate our uncertainty in I]h2 before deprojection to be 
la ^ 15 Mq pc^^ (equivalent to I]h2 ~ 11-5 M© pc"^ 
after correcting for the inclination angle of the SMC). 

The following points are important to keep in mind, 
since they represent limitations of our H2 map and ul- 
timately of our analysis. First, to remove spurious H2 
emission from our map we only retain contiguous regions 
of area > 4 square arcminutes (this cleans up a few islets 
of emission) and Sh2 > 23 Mq pc~^ {2a deprojected), 
and set the rest of the map to Eh2 = 0. Because of this 
and the method used to derive a local ^gdr, we cannot 
recover a pervasive H2 component. Other observations 
suggest that a large molecular component of this form i s 
very unlikely (jDickev et al.ll2000HTumlinson et al.ll2002[ ). 
thus we do not view this as a significant uncertainty al- 
though it remains a possibility. 

Second, the SMC has a complex geometry (see dis- 
cussion in i jl.ip . As a result certain lines of sight may 
contain a combination of regions with different (5gdr- 
The presence of diffuse, high-t^cDR emission along the 
line of sight will invalidate Equation [TJ The easiest way 
to correct for this is to subtract a "dust-free" component 
from the H I map to adjust the zero point of the dust- 
gas correlation. The likely magnitude of this dust-free 
H I component for the SMC is Ehi --20-40 M© pc'^ 
(|Leroy et al.ll2011| ) . We do not remove a "dust- free" com- 
ponent in the analysis presented here. Removing it would 
drive our H2 map toward lower values of I]h2- 

Third, despite our attempt at minimizing the effect 
of systematic Jgdr changes through a local determina- 
tion, variations in the dust emissivity and Sgdr between 
the dense and diffuse ISM are possible and largely un- 
constrained. If they have the expected sense of higher 
optical depths per unit gas in H2 than in H I, they will 
drive the I]h2 toward values that are too high. Thus the 
errors introduced have the same sign as those discussed 
in the previous point. Our best estimate of their com- 
bined effect is a factor of two systematic uncertaintly in 
Eh2, with lower Sh2 values more likely. 

Fourth, at 98" the resolution of the H I map is some- 
what lower than that of the IR maps. Therefore at spa- 
tial resolutions better than r ~ 30 pc we have assumed 
that Shi is s mooth at scales below th e 1.6' resolution of 
our H I map (|Stanimirovic et aDllQQQf ) . This assumption 
is almost certainly flawed in detail, but will average out 
across the whole galaxy, and it has no impact on quan- 
tities determined on r > 30 pc scales. 

Our estimate of a factor of 2 systematic uncertainty 
in our H2 determination, together with source geometry 
uncertainties discussed in i j 1.1 1 and ^3.3\ result in a com- 
bined systematic uncertainty of ~ 0.5 dex (a factor of 3) 



in the molecular depletion time in the SMC ( 



2.2. Estimating H2 from CO 

The standard practice is to estimate the amount of 
molecular gas present in a system using ^^CO observa- 
tions. This requires the use of a proportionality factor to 
convert intensity into column density or mass. We use 
the following equations 



iV(H2)=XcoI 



CO 



M„ 



"CO -^co, 



(2) 
(3) 



where our adopted proportionality constants appropriate 
for Galactic gas are Xqo = 2 x 10^° cm~^(K km s~^)~^ 
and aco = 4.4 M© (K km s~^ pc^)~-^, Ico is the in- 
tegrated intensity of the ^^CO J = 1 — ;> transition 
(usually in K km s~^), and Leo is the luminosity of the 
same transition (in K km s^^ pc^). The resulting col- 
umn density and mass are in cm^^ and M© respectively, 
and the molecular mass Mmoi corresponds to the mass 
of H2 corrected by the contribution of the cosmic abun- 
dance of He. Note that although the Galactic values are 
approximately appropria te on the small spati al scales in 
the CO-bright material (jBolatto et a l." 2008"). the y are 
most likely inappropriate on the large scales (Rubio et al.l 
19931: 1!sraei|[l997l : llsrael et al.|[200a iLerov et al.l l20Q7al 



20111) . This is not surprising: if CO and II2 are not per- 



fectly co-located their ratio will depend on the regions 
over which we are averaging. In particular, if CO is con- 
fined to the highly shielded high surface density cores 
while H2 is more w idespread fe.g.. iLequeux et al1ll994l : 
iBolattoet a"Llll999D . Xco will be small on CO-bright re- 
gions and large when measured on large spatial scales 
(CO freeze-out into grains is not important on the spa- 
tial scales considered). We discuss this further in iJH 

2.3. Tracing Recent Star Formation 

We use Ha emission, locally corrected for extinc- 
tion, to trace the surface density of recent star for- 
mation. The Ha observations we use were obtained, 
calibrated, continuum subtracted, and mosaiced by 
the Magellanic Cloud Emission L ine Survey (MCELS, 
iSmith fc The MCELS Tea^ [19991 ). Particularly for 
bright stars, the continuum subtraction may leave no- 
ticeable artifacts, which we have masked out in our anal- 
ysis. We convolve these images with the point spread 
function of the MIPS 160 /xm camera, to match them to 
the ^ 40" resolution of our molecular gas map. 

We correct the ionizing photon flux inferred from Ha 
for the effects of extinction using the MIPS 24 /ini im- 
age obtaine d by SMC-SAGE, and the prescription from 
iCalzetti et al.. (.2007.) . The implied extinctions are very 
small over most of the SMC, only becoming significant 
for the centers of the brighte st H II regions . This is con- 
sistent with the findings of I C apian et al.l ()1996[) , who 
used the B aimer decrement technique to obtain a typ- 
ical extinction Ayin ^ 0.3 mag toward b right SMC H II 
regions (see also lHarris fc: Zaritskyll200l . for a discussion 
of extinction based on fitting color-magnitude diagrams) . 
The small extinction correction should not be surprising 
since the SMC is notoriously dust poor, with a dust-to- 
gas ratio that is a factor of 5 — 10 lower than Galactic 



(jLerov et al.ll2007al) . This galaxy simply does not harbor 
much obscured star formation activity. 

The combination of Ha and 24/^m maps yields an esti- 
mate of the extinction-corrected ionizing photon surface 
density. For our analysis we convert this to a star forma- 
tion r ate surface density, Ssfr, following (jCalzetti et al.l 
120071 Equation 7), 



SSFR = 5.3 X 10-*^ pHa + O.O3IE24] 



(4) 



where Shq and S24 are respectively the surface densities 
of Ha and 24 /ini far-infrared emission, in erg s~^ kpc~^, 
and the relation to SFR assumes an underlying broken 
power-law Kroupa initial mass function. This prescrip- 
tion is local, obtained for individual star forming regions 
and calibrated against Paschen a emission. The use of 
24/xm emission to correct Ha has been tested in M 33 
by iRe laho et al. (2010) on th e scale of individual H II 
regions. iCalzetti et al.l (|2007 ^ found the metallicity de- 
pendence of the extinction correction to be ^ 20% for 
galaxies down to nebular metallicities below that of the 
SMC. 

The very high spatial resolution of our data, r ~ 12 pc 
(though much longer along the line of sight) adds some 
complication to the concept of a star formation rate. The 
calibration by Calz etti et al.l ()2007[ ) is derived on spatial 
scales of hundreds of parsecs. On the small scales cor- 
responding to our full resolution any individual line-of- 
sight may poorly sample the high mass end of the ini- 
tial mass function, be populated by stars resulting from 
a single star formation episode, or emit Ha radiation 
resulting from ionizing photons originating in an adja- 
cent region. Spatial smoothing to larger scales removes 
these concerns. Throughout this paper we consider spa- 
tial scales ranging from r ~ 12 pc to the whole SMC. At 
the smallest of these scales, the idea of local star forma- 
tion rate may break down, but at resolutions of 200 pc or 
1 kpc (and certainly for the whole SMC) the application 
of a standard Ha to SFR conversion should be adequate. 
These scales also allow a fairer comparison to most other 
extragalactic studies. 

It is important to note that in the results presented in 
the following sections we do not remove a diffuse ionized 
gas component, which may be as high as 40% of the to- 
tal Ha emission in the SMC (Kennicutt et al. 1995). The 
200 — 600 pc scales on which .Calzetti et al.. (,2007! ) per- 
formed their calibration are much larger than the ^ 12 pc 
on which we carry out our study. Thus it seems likely 
that part of the diffuse ionized emission due to the escape 
of ionizing photons from H II regions is already included 
in the calibration, suggesting that 40% is an upper bound 
to the SFR correction due to diffuse emission. To assess 
the impact of the diffuse Ha on our results we performed 
an analysis where we spatially filter the Ha image on sev- 
eral scales to rem ove the diffuse compo nent, along the 
lines described in IRahman et al.l (|2011L and references 
therein). We found that the results presented in the fol- 
lowing sections are robust to the presence of diffuse Ha. 
The main effect of removing a diffuse ionized component 
is to correspondingly lengthen the gas depletion time, 
or equivalently somewhat lower the local star formation 
efficiency. 

3. RESULTS AND DISCUSSION 



Maps of S]h2, Sgas, and Ssfr allow us to study sev- 
eral aspects of the star formation law. In the following 
sections we will characterize the global properties of the 
SMC, compare molecular gas and star formation (the 
"molecular star formation law"), total neutral gas and 
star formation (the "total gas star formation law" ) , and 
measure the molecular fraction as a function of surface 
density. 

3.1. Global Properties 

Our study allows us to derive a number of interest- 
ing integrated properties for the SMC. Integrating our 
extinction-corrected Ha map, we obtain a global star for- 
mation rate SFRsmc ^ 0.037 M0 yr~^. The extinction 
correction accounts for ^ 10% of the global SFR. About 
30% of the integrated Ha emission arises from extended 
low surface brightness regions, with equivalent star for- 
mation rate densities of Shq < 5 x 10"'^ M© yr~^ kpc~^, 
likely reflecting long mean free paths for ionizing photons 
escaping H II regions (with perhaps some contribution 
from other sources of ionization). The global SFR we 
measure is very comparable to the present SFR obtained 
from the ph otometric analysis of the r esolved stellar pop- 
ulations by iHarris fc Zaritsk-yl ()2004D , and only slightly 
lower than t he SF Rsmc ^ 0.05 M0 yr~^ obtained by 
iWilke et al.l (|200l based on the study of the SMC far- 
infrared emission. Note, however, that in the latter case 
the estimate is based on correcting up considerably the 
standard FIR calibration by assuming a much smaller 
dust optical depth in the SMC. It also relies on the syn- 
thetic model grid of Leitherer & Hcckman (199^), which 
has b een superseded by th e newer models incorporated 
in the ICalzetti et ahl (|2007D calibration. In any case, this 
range of values is representative of the typical uncertain- 
ties in determining SFRs, also applicable to our value. 

Integrating the molecular mass over the regions of the 
map with signal above our 2cr threshold of 23 M© pc~^, 
we obtain M^^^ ~ 2.2 x 10^ Mq. This is a ~ 35% re- 
vision down from previou s studies using the same tech- 
nique (|Lerov et al.ll2007a[) . mostly due to differences in 
the estimation of the dust-to-gas ratio and well within the 
factor of '^ 2 systematic un certainties presen t in the anal- 
ysis (for discussions see iLerov et al.ll2007al . [20091. The 
global CO lumi nosity of the SMC is approximately 1 x 10^ 
K km s~^ pc^ (jMizuno et al.ll200ll assuming a 20% cor- 
rection for flux in unmapped regions). This yields a glob- 
ally averaged conversion factor from CO to molecular 
mass of aco ~ 220 M© (K km s^^ pc^)^^. We return 
to aco in the next section. 

The gas depletion time, Tdop, is a convenient way to 
parameterize the normalized star formation rate: Tdcp is 
the time needed for the present SFR to exhaust the ex- 
isting gas reservoir (Young et al. 1986). Sometimes the 
inverse of this quantity is presented as a so-called star 
formation efficiency (SFE), indicating the fraction of the 
gas reservoir used in 10^ yr. Galaxies with high SFE can 
only sustain their present star formation activity for a 
short period of time, and thus are experiencing a star- 
burst. 

From this integrated SFR and Mmoi, the global molec- 
ular gas depletion time of the SMC is approximately 
T™°i = M^oi/SFR ~ 0.6 Gyr. For comparison, the 
global molecular depletion time inferred from the CO 



luminosity and a Galactic conversion factor (aco = 4.4) 
would be ~ 0.01 Gyr. Although a Galactic conversion 
factor is clearly inappropriate, this exercise highlights the 
need for an alternative tracer of H2. The CO in the SMC 
is either dramatically underluminous for the observed 
level of star formation activity, or the SMC is in the midst 
of a starburst that will exhaust the molecular gas reser- 
voir in only ~ lO'^ yr. Although the s tudy of the global 
star fo rmation history of the SMC by [Harris fc Zaritskvl 
(j2004( ) finds its present SFR to be somewhat larger than 
the past average, the magnitude of that effect cannot ex- 
plain such an implausibly short molecular gas depletion 
time. 

Given its atoniic m ass (Mhi — 4.2 x 10* M©; 
iStanimirovic et all 119991 ) . the SMC keeps only ~ 5% of 
its gas in the molecular phase. Thus it appears gen- 
uinely poor at forming GMCs. The corresponding total 
gas depletion time is t^^^ ^11.8 Gyr, so the SMC has 
enough fuel to sustain its current rate of star formation 
for approximately a Hubble time. By comparison, the 
median total gas depletion time for the 12 large galaxies 
in the THINGS/HERACLES sample is t^^^ « 6.0 Gyr 
(jLerov et al.l 120081) . Although a wide range of depletion 
times are present, 1 < t|^p < 14 Gyr, 10 of these galaxies 
have T|gp < 8 Gyr. In contrast, in H I-dominated dwarf 
galaxies and the outer disks of spirals t^^^ > 10 Gyr 
(|Bigiel et al.ll2010bl) . 

3.2. Relation Between Molecular Gas and Star 
Formation 

Figure [2] shows the correlation between molecular gas 
and star formation activity traced by the extinction- 
corrected Ha at several spatial resolutions: r ~ 12 pc 
(grayscale), r ~ 200 pc (red squares, binned), and 
r ~ 1 kpc (black circles). The gray area and the corre- 
sponding white contours show that a fairly well-defined 
relation exists between the star formation activity and 
the surface density of molecular gas obtained from the 
dust continuum modeling described in ^2.l\ even on our 
smallest scales. 

Note that the x and y axes of this plot are not implic- 
itly correlated. The abscissa contains information from 
the 70, 100, and 160 /im far-infrared continuum as well 
as the H I map, while the ordinate is chiefiy Ha. The 
small extinction correction derived from the 24 /im data 
has a negligible effect on the correlation. Also note that 
because we correct for dust temperature when deriving 
the dust surface density S„ioi should, in principle, also 
be independent of heating effects. 

The molecular gas depletion time depends on the scale 
considered (Figured)). On the smallest scales considered, 
r ^ 12 pc, the depletion time in the molecular gas is 
log[r^°V(lGyr)] ~ 0.9 ± 0.6 (t^°i - 7.5 Gyr with a fac- 
tor of 3.5 uncertainty after accounting for observed scat- 
ter and systematics involved in producing the H2 map 
as well as the geometry of the source). As mentioned 
in the previous section, t^"°' shortens when considering 
larger spatial scales due to the fact that the Ha and H2 
distributions differ in detail, but are well correlated on 
scales of hundreds of parsecs (Fig. [T|). On size scales 
of r ~ 200 pc (red squares in Figure [2]) , corresponding 
to very good resolution for most studies of galaxies be- 
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Fig. 2. — Molecular star forination law 
scale shows the binned two-dimensional dii 
to S„joi correlation at a resolution r ~ 12 pc, where ^„iol is derived 
from the dust modeling. The intensity scale is proportional to the 
number of points that fall in a bin, with white contours indicating 
levels that are 20%, 40%, 60%, and 80% of the maximum. The 
vertical gray dashed line indicates our adopted 2a- sensitivity cut 
for the r ~ 12 pc data. The red squares and black circles show the 
results after spatial smoothing to r ~ 200 pc and r ~ 1 kpc, respec- 
tively (the sensitivity limit has been moved down accordingly, as- 
suming Gaussian statistics) . The bars in the r ~ 200 pc data show 
the standard deviation after averaging in Smol bins. The black 
contours, placed at the same levels as the white cont ours, show the 
distribution of S^qj derived from CO observations I IMizuno et al.l 
120011) using the Galactic CO-to-H2 conversion factor. The dotted 
lines indicate constant molecular gas depletion times t^°^ = 0.1, 
1, and 10 Gyr. The dashed line indicates the typical depletion 
time for normal disk galaxies t^°^ ~ 2 x 10^ years IIBigiel et al.l 
120081) . The t^°' in the SMC is consistent with the range observed 
in normal disks for E^ioi derived from dust modeling. 



yond the Local Group, the molecular depletion time is 
logkdcp/(lGyr)] - 0.2 ±0.3, or t^°^ « 1.6 Gyr. The de- 
pletion time on r '^ 1 kpc scales (black circles in Figure 
[2|) , corresponding to the typical resolution of extragalac- 
tic studies, stays constant for the central regions of our 
map (where the smoothing can be accurately performed) , 
log[r^°V(lGyr)] ~ 0.2 ± 0.2. Thus our results converge 
on the scales typically probed by extragalactic studies. 
This constancy reflects the spatial scales over which Ha 
and molecular gas are well correlated. Although the pre- 
cise values differ, a very similar trend for rj^"' as a func- 



tion of spatial scale is observed in M 33 (jSchruba et al.l 
l2010r ). The further reduction of the depletion time when 
considering the entire galaxy {t^"^ « 0.6 Gyr) reflects 
the contribution from a component of extended Ha emis- 
sion, whic h is filtered out in the calibration of the SFR 
indicator IjCalzetti et al.ll2007l:[R"ahman et al]|2011[ ). The 
SMC is on the high end of the observed distribution of 
values for the fraction of diffuse Ha, bu t fractions of 
40% — 50% are common in galaxies (e.g., iHoopes et al.l 
11991 . 

Within the uncertainties, our results are not signifi- 
cantly different from the mean H2 depletion time ob- 
tained in studies of molecule-rich late-type disks on 
750 pc to 1 kpc spatial scales, where r^"' ~ 2.0 ±0.8 Gyr 



averaged over regio ns with molecular emission (SFE w 
5% within 1 Gy r, IBigiel et all [20081: iLerov et all [20081: 
iBigiel et al.ll201ll) . This is the methodology used in re- 
solved studies of r^"' in more distant galaxies, and our 
results are directly comparable. 

This contrasts sharply with the results using the CO 
map ( black contours') obt ained by the NANTEN tele- 
scope (iMizuno et al.l[200ll ). using a Galactic CO-to-H2 
conversion (the assumption in many extragalactic stud- 
ies). The CO distribution is offset by a factor of ^ 40 
from the H2 distribution. This offset corresponds to the 
most common aco implied by our dust-map, aco ~ 185 
Mq (K km s~^ pc^)~^(very similar to the global aco)- 
Both values of aco are broadly consistent with CO-to- 
H2 conversion factors obtained by previous dust contin- 
uum modeling and virial mass techniques on large scale s 
(|Rubio et al.lflQOl [!srael[l997l: ILerov et al.ll2007aL[20TT[) . 
though factors of 2-3 discrepancy persist from study-to- 
study. They differ, however, from estimates based on 
small-scale virial masses toward the CO-bright peaks, 
which tend to ob t ain values of ag g closer to Galactic 
([Israel et all [200I iBlitz et al.|[20OT: iBolatto et al] [20081: 
iMiiller et al.ll2010l) . This discrepancy between CO-to- 
H2 conversion factor on the large and on the small 
scales can be understood in terms of the existence and 
mass dominance of large molecular envelopes poor in 
CO. Such envelopes are expected at all metallicities 
([Glover fc Mac Lowll2010l) . and at the low metallicity of 
the SMC they hkely constitute the dominant reservoir of 
molecular gas ([Wolfire et al.ll2010l) . 

3.3. Relation Between Total Gas and Star Formation 

Figure [31 shows the total gas star formation law for 
the SMC, the relation between Esfr and total (H I -|- 
H2) gas surface density Sgas, as well as Ssfr vs. Ehi 
(which is almost the same, as atomic gas dominates). 
This relationship may be more complex than the molec- 
ular star formation law, resulting from a combination 
of phase balance in the ISM and the relative efficien- 
cies of different types of gas at forming stars. Using 
a power-law ordinary least-squares bisector fit we find 
that log(EsFR) = (2.2 ± 0.1) log(Egas) + (-6.5 ± 0.1) for 
Egas > 10 M0 pc~^ at the full spatial resolution of the 
observations. This is very similar to the typical 1 +p « 2 
slope measure d for this relation in H I-dominated regions 
of galaxies (Bi giel et al.ll2010bl) , or the typical power-law 
index of the Ss fr to Eht relatio n observed in faint dwarfs 
bv .Rovchowdhurv et al.l (|2009D . 

Thus the relation between Esfr and Egas is steep in the 
SMC, and similar to Esfr vs. Egas at low surface density 
in normal galaxies, but the relation in the SMC is notice- 
ably displaced toward larger total gas surface densities 
compare d to large spirals (this was already evident in the 
study bv [Kennicutt et al.|[l995[ ). The atomic surface den- 
sity in normal metallicity galaxies almost n ever exceeds 
a sat uration point of Shi ~ 10 M© pc~^ ([Bigiel et al.l 
[2008") averaged over ^ kpc scales, although smaller-scale 
higher-column H I clouds are observed in the Milky Way 
(e.g., iHeiles fc Trolandl[2003( ). That is not the case in 
the SMC, where Ehi reaches values of ~ 100 Mq pc~^ 
(A^(II I) ■^ 1 X 10^^ cm~^, see the white contours in Figure 
[3]). This is not purely the result of the high spatial reso- 
lution; high Ehi persists even averaged over large spatial 



scales ([Stanimirovic et aT1ll999( ). It may be partially en- 
hanced by the complex geometry of the source; the SMC 
is an interacting galaxy that may have significant elon- 
gation along the line of sight , not a simple flat disk (e.g., 
lYoshizawa fc No guchi 2003") . All surface densities (av- 
eraged over large scales) would be reduced by a factor 
cos i/0.77 if the geometry is disk-like and inclination ex- 
ceeds 40°, or possibly by a larger factor if the galaxy 
is also elongated along the line of sight. These correc- 
tions are not large enough to explain the full magnitude 
of the effect we see, in light of t he H I kinematic analy- 
sis by [Stanimirovic et aJT ([2004D , but could conceivably 
contribute a factor of ~ 1.5 — 2 (where the 50% factor 
corresponds to changing the inclination from i = 40° to 
i = 60°, and a further 20%-30% factor is an estimate for 
elongation, obtained by evaluating the contribution from 
gas at high velocities to Ehi, see §3.4|) . 

This displacement of appreciable star formation activ- 
ity towards high gas surface densities means that ap- 
plying the Kennicutt-Schmidt relation observed for nor- 
mal galaxies to the total gas surface density in the SMC 
would lead to a dramatic overprediction of the star for- 
mation rate. For example, in a normal disk we would 
expect total gas surface densities of Egas ~ 20 Mq pc~^ 
to be associated with star formation r ates of Esfr. ^ 0-01 
Mq y r"^ kpc-2 on ~ 1 kpc scales ([Bigiel et al.[ [20081 
[201 li) . while in the SMC the typical SFR associated with 
such gas surface density is one to two orders of magnitude 
lower. Conversely, using the observed SFR to estimate 
gas content (as it is sometimes done in studies of the high 
redshift universe) would lead to dramatic underestimates 
of the amount of total gas present. 

This offset in the total gas star formation law stands 
in stark contrast to the conclusions in §3.21 about the 
relation between star formation activity and molecular 
gas, where the SMC appeared very similar to normal 
disk galaxies. The fact that the molecular SFE in the 
SMC is within the range observed in other galaxies im- 
plies that, within the usual (factor of 3) uncertainties, 
the observed star formation rate accurately reflects the 
amount of molecular gas present. 

This difference implies that there are more quantita- 
tive similarities between the distributions of molecular 
gas and recent star formation than between atomic gas 
and recent star formation. We plot this directly in Fig- 
ure [31 which compares the cumulative distributions of 
star formation and gas surface density at ~ 200 pc reso- 
lution, a few times larger than the size of a large GMC 
in the Milky Way. This represents a typical size scale 
over which gas and star formation should be correlated. 
The abcissa in Figure [H corresponds to the fraction of 
the total SFR below a particular value of the Esfr- The 
ordinate corresponds to the fraction of gas accumulated 
in those lines of sight. 

Figure 131 shows that star formation and molecular gas 
track linearly with each other, while the nonlinear shape 
of the Ehi distribution reflects the fact that most of 
the H I is found on lines of sight with little star for- 
mation activity. A similar phenomenon is observed in 
other faint dwarf galaxies ([Roychowdhurv et al.l [2009L 
[20111) . This lack of correspondence between H I and 
SsFR is even more accute at r '-^ 12 pc, while H2 and 
star formation continue to track each other well. Ap- 
proximately 85% by mass of the H I in the SMC is in 
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Fig. 3. — Total gas star formation law in the SMC. The gray scale 
shows the two dimensional distribution of the correlation between 
SsFR and Sgas, where Sgas is the surface density of atomic plus 
molecular gas corrected by Helium. The white contours indicate 
the correlation due to atomic gas alone, which dominates the gas 
mass (and Sgas) in the SMC. The contour levels, and the dotted 
lines indicating constant t^^"^, are at the same values as in Figure 
[2] The d ash-dotted, dashed, and solid lines indicate the loci of the 
model bv lKrumholz et al.l Il2009d ') (KMT09) for dumping-factor by 
mctallicity products cZ = 5, 1, and 0.2 respectively. The first two 
bracket the behavior of most galaxies observed at 750 pc resolu- 
tion (see KMT09, Figure 1), while cZ = 0.2 would be the value 
expected for the SMC with unity dumping-factor (a reasonable 
assumption for the spatial resolution of the observations presented 
here, r ~ 12 pc). Note that the surface density at which H I starts 
to saturate in the SMC is Ehi ~ 50 Mq pc~^ (the typical sur- 
face density is Shi ~ 85 Mq pc"'^ at the full resolution of the H I 
data), considerably larger that the typical value in normal metal- 
licity galaxies where Shi < 10 Mq pc'^ l|Bigiel et al.ll2008l . [20Til) . 
As a consequence any use of a "standard" total gas-star formation 
correlation for the SMC would dramatically underpredict total gas 
surface densities, or overpredict star formation activity. This is not 
true for molecular gas, as we discuss in the previous figure. 

the warm phase, according to observation and model- 
ing of 28 hnes-of-sight tow ard background continuum 
sources (Dicke y et aLl I2000D . Together with Figure |31 
this strongly suggests that most H I in the SMC is not 
directly related to the star forming gas, i.e., it does not 
belong to atomic envelopes of molecular clouds. 

3.4. Molecular Gas Fraction 

Figure [5] shows our final basic observational result, 
the molecular fraction of the SMC as a function of to- 
tal gas surface density. Grayscale shows the density 
of our data in (I]H2/SHi)-(Sgas) space at 200 pc res- 
olution, the blue contours show it at 12 pc resolution. 
The curved labelled lin es indicate model predictions by 
iKrumholz et al.l ()2009d) . At ^ kpc resolution, most mas- 
sive star-forming disk galaxies lie between the curves la- 
belled cZ = 1 and cZ = 5 (see ^4.ip . For a given total gas 
surface density, the SMC has molecular fractions much 
lower than these large galaxies with the offset often an 
order of magnitude or more. We discussed in §1. li the fact 
that the complex distribution of H I in the SMC along 
the line-of-sight is a source of uncertainty. We can obtain 
a rough estimate of the effects of H I not associated with 
the disk of the SMC on the molecular ratio by recalculat- 
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Fig. 4. — Cumulative distribution function of gas mass in lines-of- 
sight with increasing star formation activity at a spatial resolution 
of 200 pc, plotted against the cumulative distribution function of 
SFR for the same lines-of-sight. The abcissa corresponds to the 
fraction of the total SFR below a particular value of the Ssfr. The 
ordinate corresponds to the fraction of gas (H I in black, H2 in gray) 
accumulated in those lines of sight. The squares show the locations 
of particular values for log[SspH,], in units of Mq yr~^ kpc~^. For 
example, about 28% of all the extinction-corrected Ha emission 
in the SMC is found in lines-of-sight with corresponding surface 
densities Ssfr < 10~^'^ Mq yT~^ kpc~^, and those lines-of-sight 
contain ~ 69% and ~ 17% of all the atomic and molecular mass, 
respectively. This plot highlights the fact that most H I is found 
in regions with little star formation, even on scales of 200 pc. 

ing Shi after removing H I emission outside the velocity 
range visr ~ 120 — 180 km s~^, taken to be representa- 
tive of the disk of the galaxy. We find this exercise lowers 
Ehi by at most 30% the faint regions of the Wing, and 
more typically ~ 10% — 15%. This is a small correction 
in the molecular ratio, well within our uncertainties in 
Eh2 alone, and although it may significantly contribute 
to the observed dispersion in Eh2/Ehi it cannot be the 
cause of its offset with respect to normal galaxies. The 
SMC is strikingly bad at turning its wealth of atomic 
gas into molecular gas, particularly given the very high 
surface densities found in this galaxy. 

3.5. Synthesis of Results 

Using our dust-based Eh2 rnap, we showed that to first 
order the molecular star formation law in the SMC re- 
sembles that in disk galaxies. There is still room within 
the uncertainties for a factor of 2-3 decrease in r^p , but 
our best estimates at 0.2-1 kpc resolution imply very 
good agreement between this low-metallicity dwarf and 
more massive disk galaxies. Note that since the scaling 
is linear, this result is insensitive to uncertainty in in- 
clination or other aspects of the SMC's geometry. By 
contrast, the total gas star formation law is offset signif- 
icantly from that observed in large galaxies. The SMC 
harbors unusually high Ehi and low Esfr at a fixed total 
gas surface density (although the Esfr vs. Egas distri- 
bution moves closer to the loci of large spirals if the star 
and gas are in a disk inclined by i > 40° , or if the galaxy 
is elongated along the line of site). At a given Egas, the 
molecular gas fraction is also offset to values lower than 
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Fig. 5. — Ratio of molecular to atomic gas in the SMC. The 
blue contours and the gray scale show the two-dimensional distri- 
bution of the ratio Sh2/Shi versus Sgas on scales of r ~ 12 pc and 
r ~ 200 pc respectively (note that the hard edge present in the blue 
contours at low ratios and low Sgas is the result of our adopted 2cr 
cut in Sh2)- The dotted horizontal line indicates Sh2/Shi = 1, 
denoting the transition between the regimes dominated by H I and 
H2. The dash-dotted, dashed, and solid lines show the predictions 
of KMT09 for different values of the cZ parameter, as in Figure 
[31 For the SMC, cZ = 0.2 at r ~ 12 pc, and the KMT09 curve 
overestimates the molecular to atomic ratio by a factor of 2 — 3. 
At r ~ 200 pc we may expect cZ Ri 1 using the standard clumping 
factor c = 5 adopted by KMT09 for unresolved complexes. Al- 
though molecular gas in the SMC is highly clumped, the atomic 
gas is not, so the cZ = 1 curve overestimates Sh2/Shi at 200 pc. 
The thick gray and black contours indicate the predicted surface 
density ratio of gas in gravitationally bound complexes to diffuse 
gas, Egbc/S(jiff , in OMLIO and in the model modified to incorpo- 
rate the extra heating of in the diffuse gas (OMLlOh, see ^4.2.11 1. 
respectively. The contours are calculated for the metallicity and 
distribution of stellar plus dark matter density in the SMC. The 
original OMLIO prediction for Sg^c/Sjiff is considerably higher 
than the observed Sh2/Shi. The tightness of the contours is due 
to the fact that the self-gravity of the gas dominates over the stellar 
plus dark matter contribution, thus there is an almost one-to-one 
correspondence between Sgas and the prediction for Sg^c/Sjiff . 



those observed in massive disk galaxies, by typically an 
order of magnitude. 

Two natural corollaries emerge from these observa- 
tions. First, molecular clouds in the SMC are not ex- 
traordinarily efficient at turning gas into stars; star for- 
mation proceeds in them at a pace similar to that in 
GMCs belonging to normal disk galaxies. This sug- 
gests that, down to at least the metallicity of the SMC 
[Z ~ Zo/5), the lower abundance of heavy elements does 
not have a dramatic impact on the microphysics of the 
star formation process, although it does appear to have 
an important effect at determining the fraction of the 
ISM capable of forming stars. 

This is not a foregone conclusion. For example, it 
is conceivable that the low abundance of carbon and 
the consequent low dust-to-gas ratio and low extinction 
would affect the ionization fraction in the molecular gas. 
This may result in changes in the coupling with the mag- 
netic field, perhaps slowing the CMC collapse and result- 
ing in lower SFE and longer t™^. Or, alternatively, the 

low abundance of CO (an important gas coolant in dense 



molecular cores) could make it difficult for cores to shed 
the energy of gravitational contraction, slowing their col- 
lapse and again resulting in lon ger time scales for con - 
suming the molecular gas (but see lKrumholz et 3111201 It ) . 
Our result implies that to first order metallicity does not 
have a large impact on the rate at which star formation 
proceeds locally in molecular gas in the SMC. Firming 
up this conclusion, however, requires detail ed studies of 
molecular cloud lifetimes (for example, see iFukui et"al] 

[Mi). 

Second, these observations provide very strong evi- 
dence that star formation activity is related directly to 
the amount molecular gas, with H I coupled to SFR only 
indirectly. This should be tempered by the considera- 
tion that, as pointed out in i )2.11 our dust-derived H2 
map may include a contribution from very cold, strongly 
self-absorbed H I (such as that sometimes associated with 
molecular cloud envelopes) which we cannot easily dis- 
entangle from molecular material in our analysis. The 
strong relation between molecular material and star for- 
mation explains some puzzling o bservations in t h e con - 
text of H I dominated systems. iWolfe fc ChenI (|2006f l 
searched for low surface brightness galaxies in the Hub- 
ble Ultra Deep Field, with the expectation that, based 
on Damped Lyman Alpha (DLA) neutral gas column 
density statistics and the star formation law, a mea- 
surable fraction of the sky should be covered by low 
brightness objects if star formation takes place in DLAs. 
Their study shows that star formation activity in DLAs 
is suppressed by over an order of magnitude (factors of 
30 — 100) with respec t to the predicti ons based on the to- 
tal gas Schmidt law. IWolfe fc ChenI indeed suggest that 
part of the explanation m ay be a low molecular fr action 
in DLAs. More recently, iKru mholz et al.l (|2009a[ ) show 
that although the observed distributon of column densi- 
ties in QSO-DLA systems theoretically requires the ex- 
istence of a significant cold phase, they are inconsistent 
with the expectation of large molecular fractions. This 
is simply a reflection of the fact that, given their low 
metallicities, their densities are not large enough to sus- 
tain significant molecular fractions. The analysis of the 
SMC presented here shows that star formation activity 
is directly proportional to the molecular content. 

4. COMPARISON TO STAR FORMATION MODELS 

How do our results in the SMC compare with pre- 
dictions from models? The last few years have seen a 
range of very important theoretical and computational 
modeling effor t concerning the star formation law in 
galaxies (e.g.. iSchavd 12004 iBlitz fc RosolowsJwl 120061: 



Schaye & Dalla Vecchia' '2008^ [Robertson fc Kravtsovl 
2008; Gncdin & Kravtsov 2010130, to mention just a 
few) . In the following sections we will focus on two recent 
theoretical models for the dependence of star formation 
and phase balance o n local gas content and o ther galactic 
properties, t hose of IKrum holz et al.l ()2009c f). henceforth 
KMT09, and lOstrikeret al. (201^ henceforth OMLIO. 
Since the models assume simple geometries, the reader 
should keep in mind the caveats raised in §1.1! about the 
complex line-of-sight geometry and overall structure of 
the SMC throughout this section. 

The models of KMT09 (summarizing a series of pa- 
pers) and OMLIO adopt different simplifications and fo- 
cus on different aspects of gas phase balance in the ISM, 
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and are therefore complementary. The KMT09 model 
adopts the simplification that all neutral gas is gathered 
into high-density, cold atomic-molecular complexes with 
surface density Scomp a factor c larger than the mean gas 
surface density Egas averaged over large scales (radio ob- 
servations typically sample scales > 0.3 kpc for galaxies 
outside the Local Group). The parameter c is not pre- 
dicted by KMT09, but the comparison with observations 
suggests c ^ 5 on scales r ~ 750 pc, and by definition 
c '^ 1 on scales approaching the size of molecular cloud 
complexes where Ecomp dominates the surface density of 
gas (several tens of parsecs to ^ 100 pc for the Milky 
Way). Warm, diffuse H I gas is assumed to represent a 
negligible fraction of the total gas content. This assump- 
tion does not in fact appear to be satisfied globally in the 
SMC, based on our estimated H2 mass of 2.2 x 10^ M© 
and an H I mass of 4.2x 10^ Mf^, most of which is believed 
to be warm (jDickev et a l. 2000). The cold H I-f H2 com- 
ponent, however, probably still dominates over warm H I 
locally. In the KMT09 model, the H2/H I balance within 
complexes is determined by shielding of dissociating radi- 
ation, and depends primarily on ZS, 
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Star formation 
is assumed to take place exclusively within gas which is 
H2 rather than H I, at the typical rates and efficiencies 
observed locally, of a few percent per free-fall timescale. 
Comparison between our observations and KMT09 are 
better carried out at the full resolution of the data, at 
which the dumping-factor c should approach unity. 

The OMLIO model adopts the simplification that all 
neutral gas is divided between a diffuse component (con- 
sisting of both warm and cold atomic gas in mutual 
pressure equilibrium), and a gravitationally-bound com- 
ponent (consisting of cold atomic and molecular gas in 
unspecified proportion). The amount of diffuse atomic 
gas is set by the requirement that heating (primarily 
FUV) balances cooling, with the mean density of the 
diffuse medium (and hence the cooling rate) set by dy- 
namical equilibrium in the total gravitational potential 
(provided by stars, gas, and dark matter). Star forma- 
tion, which produces the FUV, is assumed to occur at 
a constant rate within the gravitationally-bound com- 
ponent. Because only the very densest, highest-column 
gas within a gravitationally-bound cloud actually forms 
stars and the distribution of densities and columns in 
highly-turbulent clouds depends on temperature but not 
H2 content, OMLIO assumed that the star formation rate 
in the gravitationally-bound component is independent 
of the large-scale H2/H I proportions within this compo- 
nent, with a depletion time of 2 Gyr. Comparison to ob- 
servations should be carried out on scales where there is a 
mix of phases in equilibrium within a resolution element, 
which we take to be approximately correct at r > 200 pc. 

4.1. Comparison to KMT09 

The KMT09 model is very successful at reproducing 
the composite star-forming properties of local samples of 
galaxies with only three inputs: surface density of gas, 

/Egas that 
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metallicity, and a clumping factor 

approaches unity at high spatial resolution. The first two 
inputs are directly observable. On the other hand, the 
clumping factor is poorly constrained in unresolved ob- 
servations, and is unlikely to be constant if star-forming 
cloud complexes are gravitationally-bound entities that 
are isolated from their environments. The clumping fac- 
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Fig. 6. — Star formation law in the SMC at different spatial 
resolutions (compare to the full-resolution results in Figure[3)l. The 
gray scale shows the two-dimensional distribution of the correlation 
between SsFR and Egas, as in FigurefS] but now at a resolution r ~ 
200 pc. The white circles show the same information at r ~ 1 kpc 
resolution. The results of the KMT09 models for different values 
of the cZ parameter are indicated by the dash-dotted, dashed, and 
solid lines, as in Figure [3] The dotted lines indicating constant 
SFE are as in Figures [2] and [3] Degrading the resolution tends to 
move the star formation-total gas relation along a line parallel to 
the original distribution in Fig. [S] rather than shifting to the left 
following the predictions of KMT09 for different clumping factors. 



tor is introduced to correct the surface densities observed 
on large spatial scales to the "true" gas surface density 



on GMC scales, Ecomp- Fortunately for us, on the scale 
of our SMC observations the clumping factor should be 
approximately unity. Furthermore, these data make it 
possible to test the effects of resolution on measurement 
of the star formation law in H I -dominated galaxies. 

Figure [3] shows the predictions from KMT09 over- 
laid on the SMC data. We show the model results for 
three values of the clumping factor-metallicity product, 
cZ = 0.2, 1, and 5. The former value corresponds to 
the metallicity of the SMC and a clumping factor of 
unity, which is to be expected given our spatial reso- 
lution of r ~ 12 pc. The latter two values bracket 
most of the observations originally used to test the model 
(jKrumholz et al.ll2009bl ). Although the model drops off 
more steeply than the data at low Egas and overshoots 
the observations at high Egas, the overall agreement be- 
tween the SMC observations and the cZ = 0.2 model 
curve is very reasonable. 

KMT09 matches the star formation law observations 
reasonably well in most part because of its success at re- 
producing the observed molecular fraction as a function 
of column density. Figure [5] shows the distribution of the 
ratio of molecular to atomic gas as a function of total gas 
surface density at scales of r ^ 200 pc, together with the 
model tracks. Although most lines-of-sight are atomic- 
dominated, the model with cZ = 0.2 does a reasonable 
job at describing the observations. There is almost no 
differentiation in the location of different subregions of 
the SMC in this plot. Most notably, the SW end of the 
SMC Bar with the largest surface densities (where most 
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Fig. 7. — Molecular fraction Smoi/Sgas in the SMC, compari- 
son of KMT09 and measurements at 12 pc resolution. The top 
panel shows the results of KMT09 corresponding to our Sgas mea- 
surements, computed for the metallicity of the SMC and unity 
clumping factor. The bottom panel shows the measurements ob- 
tained from the H2 map presented here (Figure [TJ and the H I 
observations by IStanimirovic et al.l II1999I ) . Although the agree- 
ment between model and observations is reasonable, KMT09 tends 
to overpredict (by a factor of ~ 2) the molecular fraction (and the 
SFR) at high Sgas- This was already present in Figures |3] and [5] 
but it is much more apparent on a linear scale such as that used 
here. Note that, because of the spatial filtering properties of the 
algorithm used to produce the H2 map, our measurements are not 
sensitive to an extended low level molecular component. 



star formation takes place) tends to lie preferentially be- 
low and to the right of the model, while the NE end of 
the bar produces a lot of the scatter toward lower col- 
umn densities. The fact that KMT09 obtains reasonable 
SFRs as a function of surface density suggests that one 
of the fundamental assumptions underlying the model, 
that once gas turns molecular it forms stars with ap- 
proximately fixed efficiency and rate, is approximately 
correct. 

One important check on KMT09 is whether the clump- 
ing factor parameter will explain the SMC data at dif- 
ferent resolutions. Most of the molecule-dominated disks 
used in previous comparisons with the model have metal- 
licities that are approximately solar and are consistent 
with tracks having cZ '^ 1 — 5, suggesting clumping fac- 
tor values c ~ 1 — 5 at resolutions of several hundred 
parsecs. The results of spatially smoothing the SMC 
data for Esfr and Sgas to resolutions r ~ 200 pc and 
r ^ 1 kpc are shown in Figure [51 The gray scale and the 
white circles represent the results at 200 pc and 1 kpc 
respectively. We can see that instead of moving from the 



cZ ~ 0.2 track at r ^ 12 pc to the cZ ^ 1 track at 
r ^ 1 kpc the points tend to preserve the distribution 
of the original data, sliding down along a line with slope 
1 + p « 2 (a fit to the 200 pc resolution points finds 
1 + p — 2.1 ± 0.2). Thus even at r ~ 1 kpc the total gas 
star formation law points are found more-or-less along 
the cZ = 0.2 track. 

Why is the inferred dumping-factor parameter not 
changing much with resolution? The main reason is that 
the distribution of H I surface density is very uniform in 
the SMC. Thus going from r ^^ 12 pc to r ~ 1 kpc does 
not dramatically affect the surface density in its central 
regions. By contrast, the squares and circles in Figure 
[2] show the effect of going to larger spatial scales on the 
H2 distribution. The median of the logarithm of Eh2 on 
1 kpc scales decreases by 0.5 — 0.6 dex from its value 
at r '^-^ 12 pc, about the x-axis separation between the 
cZ — 0.2 and the cZ = 1 model tracks in Figure [6] In 
other words, if the distribution of H I were as clumpy as 
the H2, the change in spatial resolution would shift the 
points by about the separation between the model tracks 
along the x-axis. Thus, the H I and H2 gas does not fol- 
low a similarly clumped spatial distribution. This calls 
into question the assumption by KMT09 that the atomic 
medium is primarily found in shielding envelopes of im- 
phcitly cold, dense H I g as surroundin g the H 2 gas. In- 
deed, the observations of iDickev et al.l (|2000f) suggested 
that the H I in the SMC is primarily warm, diffuse gas 
rather than cold, dense clouds. 

Note that c = Scomp/Sgas, so using a constant c over 
a range of surface densities is strictly incorrect, and this 
may be partially the cause of the apparent slant of the 
data with respect to the cZ = 0.2 track, which is present 
at all resolutions. Impossing c ~ (Egas)"^, however, dra- 
matically overcorrects this effect. Fitting the observa- 
tions requires a softer correction, c ~ (Sgas)""'^, sug- 
gesting a connection between the surface density of cloud 
complexes and the density of the surrounding gas on 
large scales. It is not immediately clear why this should 
be the case in the context of KMT09. 

It is worthwhile to note that although the KMT09 
cZ = 0.2 track is in reasonable overall agreement with 
the data in Figure [SJ there is a noticeable bias toward 
overpredicting molecular ratios by factors ~ 2 — 3. Since 
the discrepancy is of the order of our claimed system- 
atic uncertainty for the H2 map, this is most meaningful 
in the sense of the relative comparison to OMLIO. Fig- 
ure [7] shows the molecular gas fraction Emoi/Sgas result- 
ing from applying KMT09 to the observed Egas distri- 
bution at full resolution (where we know the dumping- 
factor should be c ~ 1 and cZ ~ 0.2), compared to our 
measurements of the same quantity. The discrepancy 
between model and observations in the molecular frac- 
tion is most apparent for the SW tip of the Bar, which 
harbors the largest atomic and total gas surface densi- 
ties. Note that the discrepancy between the data and the 
cZ = 0.2 predictions persists although somewhat dimin- 
ished at r ~ 200 pc (Figure [5]), but the model applicable 
on those scales should have c > 1 {cZ > 0.2) since the 
molecular complexes are likely unresolved and Sgas is 
smaller than Ecomp- 

4.2. Comparison to OMLIO 
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The model of OMLIO was motivated in part by obser- 
vations indicating that the star formation rate and phase 
balance depend not just on the gas surface density, but 
also on the density of the stellar component, with the 
molecular content increasing r oughly linearly w ith the es- 
timated pressure in the ISM (jBlitz fc Rosolowskv 2006: 
iLerov et al.l 120081 ) . OMLIO explained this empirical re- 
sult in terms of the simultaneous need to satisfy ther- 
mal and dynamical equilibrium in the volume-filling dif- 
fuse component. The cooling rate in the diffuse medium 
is proportional to the gas density, which is controlled 
by vertical dynamical equilibrium (and hence can be af- 
fected by the stellar density if it is large) . The heating 
rate, associated with radiation from high-mass stars in 
galaxies with active star formation, is proportional to 
the amount of gas forming stars in gravitationally bound 
cloud complexes (GBCs). In order to balance cooling 
and heating, an appropriate partition of gas into diffuse 
and self-gravitating components is needed. The solu- 
tions of the simultaneous thermal/dynamical equilibrium 
equations were shown by OMLIO to agree very well with 
the radial profiles of star formation in a sample of spi- 
ral galaxies. The model of OMLIO docs not, however, 
make a prediction for the atomic/molecular balance in 
the ISM; in t j4.2.2[ we describe an extension of the model 
that provides this estimate, for comparison to the obser- 
vations in the SMC. 

We use the following implementation of OMLIO, based 
on their Equations 10 and 15 — 17, 



^diff — 2;Zjgas, 

^gbc ^^ (.1 ^ 2;JZjgas, 

(1 - a;)Egas 



SsFR = 



4I]SFR 

y = - X 



Zd/Zg 



Sh; 



diff - 



SSFR,0 l + 3.1(SgasVSo)0-365' 

9.5Mo pc"2 Q,y 
O.llSgbc + [0.011(Sgbc)2 + ay + lOOa/wPs, 



(5) 
(6) 

(7) 
(8) 

off) 



where x represents the fraction of gas in the diffuse 
phase, Sgas, Sdig, and Sgbc are the large-scale aver- 
ages of the surface densities of total gas, diffuse phase, 
and gravitationally-bound cold phase respectively (all in 
Mq pc~^), and psd is the midplane volume density of 
stars plus dark matter (in M© pc~'^). The parameters 
So « 10 Mq pc-2, SsFR.o « 2.5 X 10-9 Mo pc'^ yr'i, 
a ~ 5, /w w 0.5, and Tdcp « 2 x 10^ yr are respectively 
the surface density of gas and the star formation rate at 
the Solar circle, the ratio of total to thermal pressure, the 
fraction of diffuse gas in the warm phase, and the gas de- 
pletion time in the gravitationally-bound component (as 
inferred empirically from spiral galaxies). The parameter 
y is simply the normalized thermal pressure in the model, 
y = (Pth/fc)/3000K cm^"^ where k is Boltzmann's con- 
stant. For a two- phase ISM, Pth is p roportional to the 
UV heating rate (jWolfire et al.ll200"3J ). which is propor- 
tional to SsFR (as expressed in EquationUl). Finally, the 
dust-to-gas ratio and the gas phase metallicity, both rel- 
ative to the Milky Way Solar circle values, are indicated 
by Z d and Zg. For the SMC we assume Z^ = Zg = 0.2 
(e.g. , ILerov et al.. ,2007a) • 



Given a local Sgas and psd, we substitute Equations 
[5] to [8] into Equation [9l finding the value of x in the 
[0, 1] interval that satisfies it. The corresponding value 
of SsFR can then be readily found using Equation [71 

As inputs to the model, besides our total surface den- 
sity gas map, we employ the stellar surfa ce density de- 
rived from the SAGE-SMC 3.6 /xm images (iGordon et al. 



12011!) using the mass-to-light ratio from ILerov et al 



(200 8), and the dark matter density profile results from 
Bek ki fc Stanimirovid (|2009| 1. To deproject the stellar 
surface density we use a disk scale height of 2 kpc. Be- 
cause the thickness of the SMC is very poorly constrained 
this is simply a guess. Nonetheless, since the term as- 
sociated with psd in Equation [9] is much smaller than 
that associated with the gas terms, it turns out that 
the precise value of the scale height has little impact 
on the final result. Under these assumptions the typ- 
ical values for the bar and wing regions of the SMC 
are psd — 0.015 — 0.03 M© pc^^, with a dark matter 
contribution pdm — 0.008 M© pc^'^. As explained in 
OMLIO, for near-Solar metallicity and GBCs with cloud 
complexes surface densities Ecomp ~ 100 M© pc"^, indi- 
vidual GBCs are well-shielded so that I]h2 ^ 5]gbc and 
Shi ^ Sdiff. At low metallicity, however, GBCs may 
have substantial H I envelopes such that Eh2 < Sgbc and 
Shi > Sdiff'. We will revisit this issue in i i4.2.21 obtaining 
an estimate for the separate H I and H2 contributions to 

Sgbc- 

We find that following equations [S] to JH] with the same 
parameters as adopted by OMLIO, Sdiif is substantially 
lower than Shi in the SMC, yielding Sgbc/Sdiff much 
larger than the observed Sh2/Shi ratios. The distribu- 
tion of the model predictions computed for psd in the 
SMC are shown by the gray contour in Figure [5j This 
result is robust to our choice of Psd and the precise values 
of /w, Zd, and Tdep- In essence, total gas pressure (which 
is a Pth) is what balances self-gravity, and Pth must also 
be consistent with the balance between heating and cool- 
ing. The equilibrium value of Pth is insensitive to metal- 
licity (Equation [S]) because photoelectric heating is oc Zd 
whereas metal cooling is oc Zg . Consequently at the high 
column densities and pressures present in the SMC, the 
model requires that most of the gas is driven into the self- 
gravitating cold phase almost independent of metallicity, 
in order to generate the radiation field and the corre- 
sponding heating needed to attain pressure equilibrium 
in the diffuse gas. For Sdiii to approach the large values 
of Shi observed in the SMC would imply large values of 
the pressure, which necessitates either increases in Pth 
(and consequently the heating) , or a large turbulent fac- 
tor represented by a. Matching the observations with an 
increase in a requir es a very large value, a ~ 20, which 
is likely unrealistic (jBurkhart et al.l 12010). For the ref- 
erence value, a w 5, the predicted surface density of the 
diffuse phase in the SMC is Sdiff< 30 M© pc^^, in con- 
trast with the higher range of Shi evident in Figure [31 

It is important to realize that this point is pretty much 
independent of the details of the model. In any scheme 
where the gas is in pressure equilibrium and the geom- 
etry is not pathological, sustaining large diffuse-H I sur- 
face densities requires a correspondingly large heating 
term to counterbalance the strong cooling of warm H I 
at high pressure and density. The gas heating can be 
driven by the radiation field, or else dynamical energy 
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input (as in the example of increasing the a parameter) 
can puff up the disk to reduce its density and coohng 
rate. In either case, energy and momentum input from 
star formation needs to be present to provide the support 
necessary to prevent collapse. To match the typical H I 
surface density observed in the SMC at 200 pc resolution, 
Shi"^ 60 — 70 Mq pc~^, would require a thermal pres- 
sure Pth/fc '^ 24 — 32 X 10'' K cm~^ in the diffuse phase 
(assuming a = 5), far in excess of what is observed in 
the Milky Way plane. 

Attaining such pressures requires a large heating term, 
and equilibrium demands it be balanced by an equally 
large cooling. Since the cooling in neutral gas is 
dominated by far infra red fine structure lines over a 
wide range of densities (jWolfire et al.1119951 ). we expect 
that the consequence of increased heating in the dif- 
fuse medium would be observable. Note that there 
is evidence that the [C II] 157.7 /^m transition is in- 
deed bright in th e Magellanic Clouds and other low 
metallicity dwarfs (Stacev et al.lll99lUlsrael et al.|[T996l : 



iMadden et al.lll997D . but those measurements are likely 
dominated by the dense, self-gravitating phase. Testing 
this hypothesis will likely require observations of the dif- 
fuse phase. 

Before proceeding to consider the possibility of en- 
hanced heating in more detail, we note that an increase 
in the assumed inclination angle i of the SMC's plane 
would also yield a lower prediction for Sgbc/Sdiff, and 
lower Pth. For the gas-gravity-dominated case. Equa- 
tion (O yields y ~ (Sgas/9.5AfQ pc~2)2/q; for the nor- 
malized thermal pressure, while Equations dH])-® yield 
y ^ 5]gbc/(SsFR.o''dop)- If diffuse gas dominates, com- 
bining these leads to 



■^diff 



(9.5M0 pc-2)2a- 



(10) 



Taking fiducial parameter values for I]sfr,07 tac^-, and a, 
this yields Sgbc/^^gas ~ ^gas/lOOMQ pc"^. A reduction 
of deprojected Sgas by a factor /z would lower Egtc/Sdig 
by the same factor. Thus, if the typical deprojected value 
of Egas were 20 M© pc"^ rather than 60 M© pc"^, the 
predicted Sgbc/Sdig would be a factor ^ 3 lower. To 
reduce Egbc/Sdig to the level ~ 0.1 of the molecular- 
to-atomic ratio would require an inclination correction 
factor ^ ~ 0.1 (corresponding to an almost edge-on in- 
clination i ^ 85°), however, which seems unlikely. 

4.2.1. Enhanced Keating'? 

What could be the source of the extra diffuse gas heat- 
ing in the case of the SMC? In the following we explore 
the possibility of more effective diffuse gas heating hap- 
pening at low metallicities. One potential solution is 
that the local radiation field (relative to the Solar neigh- 
borhood), Gg, is enhanced at low metallicities with re- 
spect to the corresponding case at Z = 1. A simplifi- 
cation introduced in the OMLIO model, as described in 
the paper, is to ignore local variations in the propaga- 
tion of the FUV photons, which determine Gg and the 
heating of the gas. Thus, the radiation field relative to 
the Solar circle in the Milky Way is represented in the 
heating term (Equation [5]) by the factor Esfr/5^sfr,o- 
Note that the optical depth to FUV radiation, ttjy, is 
determined by dust attenuation. For identical columns 



of gas, this is lower in a low metallicity environment 
because of the lower dust-to-gas ratio Zd, such that 
ruv w l/2Ax/Av{Nii/2 x lO^^cm-^) Z^. The escape 



probability of FUV photons, /3uv — (1 



')/^uv, 



which in the limit of large tuv is /3uv ~ I/t'uVi is conse- 
quently higher at lower Zd (for a given Nh)- To account 
for this effect we can introduce a factor of /3uv ~ 1/^d 
in Equation |8l so that the normalized thermal pressure 
becomes 



y 



yzg 



4SsFR 

SSFR,0 "" TT3.1(Sgas^d/So)0-365- 



(11) 



Interestingly, there is some evidence for an enhanced 
mean radiation field in the Bar region of the SMC. Mod- 
eling of the dust emission suggests t hat the average G'n is 
~ 4 — 5 times l arger than the local iMathis et alT ()1983f ) 
radiation field dSandstrom et al.l I2010D . An extended 
warm dust component associated with the Bar is also 
observed in the recent analysis u sing Planck data by 
IPlanck Collaboration et all (1201 la|). The me an radiation 
field in the Wing bv iSandstrom et all (|2010D appears to 
be approximately Galactic, with the exception its tip and 
the star forming regions there located, despite the large 
values of Shi there present. The modeling of the dust 
emission, however, is particularly difficult in the Wing, 
where the signal-to-noise of the observations is lowest. 

With Equation [TT] replacing Equation \8\ in the model 
the resulting Sdiff reproduces the large observed Ehi in 
the SMC using the nominal (Galactic) values for the rest 
of the parameters, including a. A comparison of the pre- 
dictions of OMLIO modified for heating (OMLlOh hence- 
forth) with our measurements of the molecular fraction 
are shown in Figures [5] and El OMLlOh does very well 
at reproducing the average trends in the data, although 
local over- or under-predictions persist at the factor of 
two level. 

The above approach is perhaps unrealistically simple, 
and it is likely to overestimate the correction to the local 
interstellar radiation field due to the diminished FUV 
extinction at low metallicities. An important fraction of 
the UV radiation escapes HII regions through very low 
extinction lines-of-sight, and that fact should be taken 
into account in the FUV propagation. In more general 
terms, the local interstellar radiation field relative to the 
Solar neighborhood, Gq, is 



^/ _ (1 XjLigas J 
Tdcp^SFR.O /O 



(12) 



where ///o is the enhancement factor for the radiation 
field relative to the Solar circle value due to propagation 
effects, and it includes the escape of radiation from star- 
forming regions, as well as the propagation through the 
diffuse gas. This factor can be written as 
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./csc.O H ° „ ' (1 — /csc,o) 

'gbc,0 

1 — £'2(Tdiff/2) Tdiff,0 



1 - i?2(Tdiff,o/2) Tdiff 
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where /esc is the fraction of the UV radiation directly es- 
caping star-forming regions, E2 is the second exponential 
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integral function, Tgbc and Tdis are the optical depth to 
the UV associated with gravitationally bound clouds and 
diffuse gas, respectively, and the subscript indicates So- 
lar circle Milky Way reference values. The second term 
in square brackets, as introduced in OMLIO, asssumes 
the diffuse gas is a uniform slab with total optical depth 
Tdiff = KuvSdiff- Note that the correction factor intro- 
duced in Equation [m is just a simplification of Equation 
(fT3|) considering only the factor between the first square 
brackets, due to the gravitationally bound phase, with 
/esc = 0. Although it is possible to reproduce the ob- 
servations of the SMC using this approach, it requires 
choosing values for a number of very poorly constrained 
parameters. 

Other poorly constrained parameters in the problem 
of how the low metallicity affects the diffuse gas heating 
are the changes in the properties of the dust grains that 
couple the radiation field to the gas. The smallest car- 
bonaceous dust grains in the ISM, associated with the 
2175 A extinction bump and the mid-infrared aromatic 
features attributed t o polycyclic aromatic hyd rocarbons 
(PAHs; 'Drain e fc L]||200lt iLi fc Draind l2002f ) . are con- 
siderably deplet ed at low metallicities and in particu- 
lar in the SMC (iMadden et al.l l200l lEngelbr acht et al.l 
120081 : ISandstrom et al.ll2010[ ). If these grains dominate 
the heating, as it appears to be the case for the LMC 
(jRubin et al.ll2009l ) , the consequently lower photoelectric 
heating efficiency could make it more difficult to heat the 
diffuse gas despite the enhanced radiation field. 

We conclude that including systematic heating effects 
in the diffuse H I gas at low metallicities, due to the 
enhanced radiation escape from the dense star-forming 
phase offers a viable model, able to approximately re- 
produce the data. We will return to whether this is the 
best solution in i i4.2.3l 

4.2.2. Is all the HI diffuse? 

An alternative explanation for the observations of the 
SMC in the context of the OMLIO model is that only 
a small fraction of the observed Ehi is diffuse gas con- 
tributing to Sdiff J and the remainder of the H I is actu- 
ally in gravitationally bound structures contributing to 
Sgbc- Physically, we would expect the H I component of 
individual GBCs to be the envelopes, and the H2 compo- 
nent the shielded interiors (note that the clumpy struc- 
ture of clouds means both the "envelope" and "interior" 
may pervade the whole cloud complex). The relation 
between H 1 and II2 within GBCs would then be deter- 
mined by the physics of photodissociation described in 
the KMT09 model. In principle it should be possible to 
test this idea by searching for the kinematical signatures 
of these gravitationally bound structures in the existing 
H I datacubes. 

As noted in tj4.H we hav e argued in tJ3 3|tha t neither 
Figure m nor the results of iDickev et al.l (|2000D support 
the scenario that most of the H I is cold gas in the en- 
velopes of GBCs (which are atomic-molecular cloud com- 
plexes in the terminology of KMT09). Nonetheless, it is 
interesting to further investigate its consequences in the 
context of the OMLIO model, just as we did for the model 
of KMT09. 

To compute the amount of H I in GBC envelopes we 
can apply the molecular fraction determination from the 
KMT09 model to the self-gravitating portion of the gas. 
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Fig. 8. — Molecular fraction Sjjjol/Sgas in the SMC, compari- 
son of OMLIO and measurements at 200 pc resolution. The top 
panel shows the results of OMLlOh, the OMLIO model modified to 
include enhance d hea ting in the diffuse phase at low metallicities 
as discussed in i|4.2.1l The bottom panel shows the measurements 
obtained from the H2 map present e d her e (Figure [TJ and the H I 
observations by Stanimirov ic et al.l II1999I ). Note that, because of 
the spatial filtering properties of the algorithm used to produce the 
H2 map, our measurements are not sensitive to an extended low 
level molecular component. 



and rely on the OMLIO model described by Eqs. [5H9] 
to determine the balance between the diffuse and self- 
gravitating phases. The equations for the resulting sur- 
face densities of H I and H2 averaged over large scales 
are 



i:m = ^e.s[x + {l-x){l- 


- ./h2)] , 


(14) 


SH2=Sgas(l-x)/H2, 




(15) 



Zg. We take 
and 100 



where /h2 — /H2(Scomp, Z) is determined through Equa- 
tion 2 in KMT09. It is a function of the surface density 
of the cloud, Ecomp, and the metallicity Z relative to the 
Milky Way, which we take to be Z = Z^ 
Scomp, to be the larger of Egbc = (1 - a;)Egas 
Mq pc~^. The latter is the approximate value for the 
surface density of resolved m olecular clouds in nearby 
galaxies (" Bolatto et al.l l2008'). Thus this model assumes 
that GBCs have at least this surface density, and lower 
surface densities are the result of beam dilution in the 
observations. The fraction of all H I found in the gravi- 
tationally bound phase, /hi'^, will then be 
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Fig. 9. — Behavior of the atomic component in the OMLlOp model described in ^4.2.21 and comparison to the diffuse component in 
OMLIO, under typical conditions for the Solar circle in the Milky Way (p^d = 0.05 Mq pc~^ and Z = 1) and the SMC (psd = 0.02 
Mp pc~ and Z = 0.2). (Left) Fraction of the total H I surface density arising from the gravitationally bound phase according to Eq. 
1161 The diffuse and the gravitationally bound phases make approximately equal contributions to Shi a^t surface densities of approximately 
160 (70) Mq pc~^ for typical conditions in the local Milky Way (SMC). (Right) Surface density of H I as a function of total gas surface 
density in OMLlOp, with the prediction of KMT09 for comparison. With the contribution of H I from the gravitationally bound phase, it 
is possible to match the typical Sfji ~ 70 Mq pc~^ observed in the SMC without an increase in the heating term. 
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Because this model extends OMLIO by incorporating 
the photodissociation model of KMT09 for comparison to 
observations of the H I abundance, we henceforth refer to 
it as OMLlOp. The behavior of the H I component in this 
model is shown in Figure [S] After accounting for the H I 
in the gravitationally bound gas through the inclusion 
of the photodissociation term from KMT09, the surface 
density of atomic gas saturates at Ehi ~ 90 M© pc^^ un- 
der the typical conditions found in the SMC instead of 
the Shi ^ 30 Mq pc^^ stemming from the diffuse com- 
ponent alone. In fact, the H I associated with the gravita- 
tionally bound gas becomes the dominant contributor to 
Ehi for surface densities Sgas ^ 70 M0 pc^^ under SMC 
conditions. Inclusion of the modifications to the heat- 
ing term discussed in §4.2.11 mostlv affect this threshold 
surface density, and increase the saturation value of Shi 
only moderately (Shi ^ 150 M0 pc~^ using Equation 
[TTl for example). 

Figure [10] shows the distribution of the molecular to 
atomic ratio Sh2/Shi in the SMC according to this 
model, compared with KMT09 and the distribution of 
the measurements at 200 pc resolution as in Figure [5l 
The similarity at Sgas> lOOM© pc~2 between OMLlOp 
and the KMT09 track corresponding to cZ = 0.2 is sim- 
ply due to the fact that the H I arising from the diffuse 
component is not dominant at these surface densities. 
The inclusion of the diffuse H I in the Shi budget is re- 
sponsible for the displacement down from the KMT09 
track by a factor of ~ 0.6. 

We conclude that including the effects of photodissoci- 
ation into the calculations to compare the predictions of 
OMLIO to the H I and II2 data is physically motivated 




Fig. 10. — Ratio of molecular to atomic surface densities in the 
SMC compared to model calculations. The thick black contours 
show the predictions of the OMLlOp model discussed in ^4.2.21 
where we include the photodissociation of H2. In this model, Shi 
in the SMC has a dominant contribution from H I in the gravi- 
tationally bound phase. The prediction of the OMLlOph model, 
which includes bot h phot odissociation and enhanced heating in the 
diffuse phase (see i|4.2.4[l , is also shown. Lines from KMT09 are as 
in Figure [5] 

and potentially offers an alternative explanation for the 
observed low values of the Sh2/Shi ratio. As we discuss 
in the next section, however, we think this is not the 
dominant consideration in explaining the data. 

4.2.3. Consequences for star formation 

We have mentioned possible tests of these ideas look- 
ing for evidence for enhanced heating or H I in the GBC 
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phase using: 1) further modehng of dust observations to 
determine the radiation field, 2) observations of diffuse 
ISM coolants such as [C II], and 3) searching for kine- 
matic signatures of gravitationally bound H I. It should 
be possible, however, to make an effort to distinguish 
between these enhanced-heating and standard heating 
variation of OMLIO model by their impact on the star 
formation law. 

In the context of OMLIO, all of the gas in the gravita- 
tionally bound component could potentially contribute 
to star formation, although due to the turbulent dynam- 
ics within GBCs, only a tiny fraction (^ 1% per free-fall 
time) is dense (and cold) enough that it actually col- 
lapses and forms stars. In the variant of the model with 
enhanced heating, OMLlOh, the gravitationally bound 
phase is depressed to allow more diffuse Shi. In the ex- 
tension of the model that includes an estimate for the 
contribution to Shi from photodissociated envelopes of 
GBCs, OMLlOp, the gravitationally-bound component is 
about as abundant as at higher metallicities, and the low 
observed Sh2/Shi ratio is explained by limited shielding 
of the gravitationally bound gas at low metallicity (as for 
the dense atomic- molecular complexes in KMT09). As a 
consequence of the higher abundance of gravitationally- 
bound gas, the SFR expected in OMLlOp at a fixed gas 
surface density is significantly higher than in OMLlOh. 

Figure [11] shows the SFR predicted by OMLIO 
with and without additional heating (taking psd = 
O.O2M0 pc-^and Z = 0.2 for the SMC), compared to 
the observations smoothed to 200 pc and 1 kpc resolu- 
tion and to the KMT09 model, as in Figure [6] Because 
the specific star formation rate in GBCs is assumed to 
be insensitive to their exact atomic/molecular balance, 
there is no difference in the predicted star formation rate 
between OMLIO and OMLlOp. Similarly, OMLlOh and 
OMLlOph (see below) have the same SgpR because they 
have the same Sgbc for a given Egas- The predictions for 
SsFR from OMLIO, OMLlOp, and KMT09 with c = 5 
are quite similar for Sgas^ 2OM0 pc^^, because for all 
of these models in this range of Sgas, the majority of 
gas is concentrated in cold, dense, star-forming atomic- 
molecular complexes with Ecomp ^ lOOM© pc~^. These 
models overpredict the Esfr observations by about a fac- 
tor of 6. On the other hand, OMLlOh (and OMLlOph) 
do a very good job at matching not only the average 
level of SFR activity but also the slope of the EspR-Egas 
relation, despite the simplicity of the heating correction. 

4.2.4. A unified model 

It is possible to think of further modifications 
to OMLIO. For example, we can extend the 

enhanced-heating model to include an estimate of the 
atomic/molecular balance within GBCs; we refer to this 
as OMLlOph. As noted above, this will produce the same 
EgpR-Egas results as the OMLlOh model as it has the 
same fraction of gravitationally bound gas, but it will 
yield a molecular to atomic ratio Eh2/Ehi a factor ~ 2.8 
times lower than OMLlOp (equivalent to a displacement 
of ^ 0.45 in the logarithm; as shown in Figure ITOl) . Such 
a model is not forbidden by the observations, and in fact 
it arguably reproduces better the lower values in the dis- 
tribution of molecular to atomic ratios than the model 
without enhanced heating. In model OMLlOph, it is also 
possible to compute the fraction of H I gas associated 
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Fig. 11. — Comparison of OMLIO predictions and total gas 
star formation law observations in the SMC. The thick solid black 
and gray lines show the results of the enhanced-heating (OMLlOh, 
OMLlOph) and original-heating (OMLIO, OMLlOp) versions of the 
model. As in Figure[6]the gray scale and white circles represent the 
SMC observations smoothed to a resolution of 200 pc and 1 kpc 
respectively. The KMT09 tracks and the constant gas depletion 
time lines are also as in Figure [6l 

with the gravitationally bound component, /hi*^! which 
has a value /^j'^ ~ 19% at a surface density Egas = 100 
M0 pc~^ , very similar to the 15% cold H I fraction mea- 
sured bv lDickevet~aLl (|2000t) . 

5. SUMMARY AND CONCLUSIONS 

We present a detailed analysis of the correlation be- 
tween the atomic and molecular gas phases in the SMC, 
and the star formation activity. This is the first time 
such analysis has been done at high spatial resolution in 
a galaxy of such a low metallicity. We us e modeling of 
the Spitzer dust continu um observations (jBolatto et al.l 
l2007tlGordon et al.ll201lD to calculate the H2 surface den- 
sities on scales of ~ 40" or r ^ 12 pc, toget her with 
observations of the H I (|S"tammirovic et al.lll999[ ) and Ha 
(iSmith fc The MCELS TeamI 119991) distributions (Fig- 
ure [T|) . An important caveat to keep in mind when con- 
sidering the interpretations of the results is the complex- 
ity of the H I distribution in the SMC, referred to in §1.11 
Another important caveat are the uncertainties involved 
in the production of the H2 map, which we discuss in 

We find that, in the regions where we detect Eh2, we 
measure a typical molecular gas depletion time t^°^ ~ 
7.5 Gyr on scales r ~ 12 pc with a factor of 3.5 un- 
certainty accounting for the scatter as well as the sys- 
tematics associated with our II2 map as well as the ge- 
ometry of the source (Figure [2|). The depletion time 
shortens when measured on larger spatial scales. On 
scales r ^ 200 pc (a typical scale over which molecular 
gas and Ha emission are well correlated in galaxies) and 
r '^ 1 kpc (a typical scale for extragalactic studies) we 
measure t^"' « 1.6 Gyr. The molecular depletion time 
for the SMC as a whole is shorter, r^°' ~ 0.6 Gyr, due 
to a large component of extended low level Ha emission. 



19 



These results are consistent with the typical depletion 
time of T^p' ~ 2 Gyr observed in normal disks on kpc 
scales (|Bigiel et all 120081 l20TTh . Consequently the rela- 
tion between molecular gas and star formation activity 
appears to be at most only weakly dependent on metallic- 
ity. This finding suggests that the molecular content can 
be used to infer the star formation activity (and vicev- 
ersa) even in galaxies that are chemically primitive and 
deficient in heavy elements. 

We also measure the relation between star formation 
and total gas surface density, which is dominated by H I 
over most of the SMC (Figure O. We find that the rela- 
tion is similarly steep in the SMC {1+p = 2.2 ±0.1) and in 
the ou ter disks of normal galaxies (l+p ~ 2; iBigiel et alj 
l2010b( ). but it is displaced toward much larger surface 
densities. In the SMC the H I surface density reaches 
values as high as Ehi ^ 100 M© pc~^, while in most 
galaxies Ehi ^ 10 M0 pc~^. As a consequence the use 
of the standard total gas to star formation rate relation 
dramatically overpredicts the star formation activity over 
much of this source. This finding supports the explana- 
tion that DLAs are deficient i n star formation beca use of 
very low molecular fractions (I Wolfe fc ChenI 120061 ) . The 
gas (or H I) depletion time for the SMC is approximately 
a Hubble time, rf^p « 11.8 Gyr. 

We compare our results with the recent analyti- 
cal models by iKrumholz et al.l ()2009d KMT09) and 
lOstriker et"all ()2010l OMLIO). For high- resolution ob- 
servations, we find that KMT09 is very successful at re- 
producing the displacement in the relation between Egas 
and EsFR toward higher surface densities (Figure IS]). We 
also find that with a clumping factor c of unity it ap- 
proximately predicts the observed relation between the 
molecular to atomic ratio Eh2/Shi and Egas (Figure [5]), 
although on a linear scale it becomes apparent that there 
is a systematic overprediction at high Egas by 0.3 — 0.5 
dex (Figure [71). The slope of Esfr vs. Egas predicted 
by KMT09 for fixed cZ is steeper than the relation ob- 
served in the SMC, both at the observed 12 pc resolu- 
tion and when the data is averaged over 200 pc and 1 
kpc scales (Figures [3] and Figure E]). Given that KMT09 
co-locates H I and H2 within the same dense cloud com- 
plexes, a strong relation between Ehi and EgpR would 
be expected at 200 pc scales, but this is not seen in the 
observations (Figure S]). We test the effect of smoothing 
the observations to larger spatial scales, and find that 
changing the clumping parameter in the KMT09 model 
does not match the effect of smoothing the observations 
(Figure [6]). The bulk of the total gas surface densities in 
the SMC stay approximately constant for spatial scales 
in the r ■^ 12 pc — 1 kpc range. This suggests that much 
of the H I is in a warm, diffuse component, which stands 
in contrast to the starting assumption of KMT09 that 
the H I consists primarily of cold gas concentrated in the 
envelopes of dense atomic-molecular complexes, with the 
warm H I a small fraction of the total ISM mass. 

In comparing the observations to the predictions of 
OMLIO, we find that the expected Egbc/Sdig is very 
large in comparison to the observed Eh2/Ehi in the 
SMC (Figure [5]), if we adopt the same heating efficiency 
and turbulence parameters for thermal/dynamical equi- 
librium as in typical large spirals with Z ^ 1. In other 
words, the model conspicuosly overpredicts the molecu- 



lar fraction in the SMC for the default parameters. This 
motivates us to investigate variations and extensions of 
the OMLIO model. In the first variation, which we call 
OMLlOh, we introduce a metallicity dependence in the 
ratio of gas heating to star formation that produces en- 
hanced heating at low metallicities. This modification is 
based on the idea that a lower dust-to-gas ratio allows 
UV to escape and travel farther from star-forming re- 
gions, potentially raising G'q significantly in the diffuse 
phase. An enhanced G'q leads to increased heating in 
the diffuse medium, permitting thermal equilibrium of 
warm H I at high pressure. Without enhanced heating, 
warm gas at the typical thermal pressure in the SMC 
{Pth/k ~ 3 X 10" K cm-3 for Ehi'- 65 Mq pc'^) would 
cool very quickly. Enhanced heating changes the bal- 
ance between the diffuse and self-gravitating phases of 
the gas. As a consequence the self-gravitating phase is 
less abundant than at normal metallicities under similar 
conditions for Egas- This model does a very reasonable 
job at matching the Eh2/Ehi derived for the SMC (Fig- 
ures [5]), although there are still local discrepancies at the 
level of ~ ±0.3 dex. 

As an extension of OMLIO, which we term OMLlOp, 
we do not change the balance between the diffuse 
and self-gravitating components, but add an accounting 
for the abundance of H I within gravitationally bound 
clouds, based on the photodissociation formalism of 
KMT09. The H I is split between the diffuse and the 
self-gravitating phases, and under the typical SMC con- 
ditions, half of the H I would be self-gravitating. The 
photodissociation calculation for the GBC component 
can also be applied to the cnhanced-heating model; we 
denote this as model OMLlOph. Model OMLlOph shares 
the reduced abundance of the self-gravitating phase with 
model OMLlOh, and includes a contribution to H I that 
originates in cold self-gravitating clouds. Note that in 
adding the photodissociation estimate for GBCs, we do 
not change the predicted SFR. Thus, model OMLIO and 
OMLlOp have the same Esfr vs. Egas as each other, as 
do model OMLlOh and OMLlOph (Figure [Il|) . 

When compared to the observed Eh2/Ehi in the SMC, 
model OMLlOp slightly overpredicts the molecular ra- 
tio (- 10 - 30% for the model vs. - 5 - 15% for the 
data; see Figure [TU)) . For the typical parameters of the 
SMC, ~ 30 - 60% of the total H I would be in GBCs for 
model OMLlOp (Figure[9]). Although this is less than in 
KMT09 (which puts all of the H I in the envelopes of cold, 
dense cloud complexes), to some extent it shares the dif- 
ficulty that in observations H I is not correlated with star 
formation activity (Figure [5]). Model OMLlOph, with a 
lower abundance of self-gravitating gas and hence Eh2 
due to enhanced heating, follows the observed Eh2/Ehi 
magnitude and slope better than model OMLlOp (Fig- 
ure fTO]) . OMLlOph also has a lower ratio of cold to warm 
H I than OM LlOp or KMT09, c onsistent with the ob- 
servations of Dickey et al.l (|2000[ ). When compared to 
star formation in the SMC, the enhanced-heating models 
(OMLlOh, OMLlOph) fit the relation between Esfr and 
Egas much better than the models that adopt the same 
heating efficiency as higher-metallicity spirals (OMLIO, 
OMLlOp), as shown in Figure [TT] 

In essence, the observed star formation rate requires 
that the eligible cold, self-gravitating gas (whether H I or 
II2) is less abundant at lower metallicities, otherwise we 
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would expect larger star formation rates for the observed 
Egas surface densities. The diminished abundance of the 
cold self-gravitating phase is supported by obse rvations 
of th e distribution of H I temperatures fPick evet alJ 
[2000), and by the relative uniformity of Shi when av- 
eraged at different scales. The possibility of an en- 
hanced radiation field that could maintain warm, diffuse 
H I at high surface density and pressure is suggested b y 
recent dust emission modeling (jSandstrom et al.ll2Q10D . 
Further observational confirmation of enhanced heating 
could come from observations of the cooling in the dif- 
fuse atomic gas. Further characterization of the fraction 
of cold atomic gas and its spatial distribution in both 
Magellanic Clouds using H I absorption/emission tech- 
niques would also be extremely valuable in constraining 
the models. 
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